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ABSTRACT 

HP9825  computer  programs  for  the  acquisition  and  analysis  of  aerosol 
data  are  listed.  The  programs  were  written  to  be  used  with  a  PMS  system 
in  wave  tank  experiments,  but  can  be  used  in  other  applications. 
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A.  INTRODUCTION 

This  report  briefly  describes,  defines  the  matrices  used  in,  and  lists 
a  series  of  programs  for  the  HP9825  used  to  accumulate  and  analyze  aerosol 
data  taken  in  conjunction  with  Particle  Measuring  Systems,  Inc.,  (PMS) 
aerosol  instruments.  The  programs  were  written  specifically  for  use  in 
wave  tank  experiments,  but  have  a  wider  application. 

The  PMS  instruments  consist  of  a  model  DAS-32  data  acquisition  system 
and  inouels  ASASP  and  CSASP  aerosol  probes  which,  between  them,  scan  the 
aerosol  spectrum  from  0.1  to  15  micrometers  radius. 

In  addition  to  a  data  acquisition  program,  programs  which  integrate 
dN/dr  spectra  over  various  radius  intervals  are  given.  The  remainder  of 
the  analysis  programs  are  used  to  examine  differences  in  spectra  taken  at 
different  times  and  to  look  at  these  differences  from  several  perspectives. 

The  programs  described  are: 

(1)  Data  acquisition 

(2)  Plots  log  (dN/dr)  vs.  log(r) 

(3)  Plots  A(log(dN/dr)  vs.  log  (r) 

(4)  Tabulates  A(log(dN/dr))  vs.  r 

(5)  Plots  log  (A(dN/dr))  vs.  log(r) 

(6)  Tabulates  1 og(A(dN/dr) )  vs.  r 

(7)  Integrates  dN/dr  over  r  from  0.2,  0.3,  0.4,  0.6,  0.8  and  1  to  7 
microns. 

(8)  Integrates  dN/dr  over  r  from  1.5,  2,  2.5,  3,  3.5  and  4  to  7 
microns. 

(9)  Plots  (l/dN1)(dN2  -  dNj),  the  fractional  change  in  dN/dr,  vs. 
log(r). 

(10)  Plots  the  fractional  change  in  dN/dr  normalized  byits  peak  value. 

(11)  Tabulates  the  aerosol  density  decay  constant  for  various  r. 
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8.  PROGRAMS 


1 .  Data  Acquisition  Program 

This  program  provides  the  interfacing  between  the  HP9825  computer 
and  the  PMS  system.  It  accumulates  data  for  a  time  period  specified  by  the 
user  and  outputs  it  to  a  printer  and  a  cassette  tape.  The  data  output  to 
the  printer  is  in  both  tabulated  and  plotted  form.  The  tabulated  data  con¬ 
sists  of  raw  counts  vs.  radius,  dV/dr  vs.  radius  and,  at  the  option  of  the 
user,  a  table  of  the  radius  bin  centers.  The  program  calculates  a  poly¬ 
nomial  fit,  whose  order  is  user  specified  (up  to  10),  of  log(dV/dr)  vs. 
log(r).  This  fit  is  plotted  along  with  the  individual  data  points  (with 
some  exclusions)  by  the  HP9871A  printer.  A  sample  output  is  shown  in 
Figure  1. 

List  of  matrices: 

Name  Function 

A[*]  Accumulates  sums  used  to  calculate  the  polynomial  fit. 

B[*]  Denotes  which  printer  symbols  to  be  used  for  the  six  (6) 
radius  ranges  in  the  plot,  except  B[7]  and  B[8]  used  to 
extrapolate  the  fit. 

C[*]  Cassette  tape  file  bookkeeping,  see  below. 

D  Current  file  number. 

E[*]  Average  radius  within  a  bin.  There  are  six  ranges  with 

15  bins  each. 

F[*]  Accumulates  sums  used  to  calculate  the  polynomial  fit. 

G[*]  Coefficients  of  the  polynomial  fit. 

H[*]  PMS  analog  housekeeping  module  outputs. 

M[*]  Used  in  formulating  extrapolations  to  polynomial  curve. 

N[*]  Radii  range  counter.  Essentially  keeps  track  of  the  total 
time  the  PMS  was  examining  a  given  radii  range. 

0[*]  dN/dr  and,  later  in  the  program,  dV/dr. 

P[*]  Sum  of  raw  counts  in  each  radius  bin  arranged  in  order  of 

increasing  radius. 
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R[*]  Defines  radii  bin  edges. 

S[*]  Sum  of  raw  counts  in  each  radius  bin  in  the  order 
specified  by  the  DAS-32. 

T [*]  dN/dr 

X[*]  Defines  which  radius  bins  in  each  range  are  to  be  excluded 
from  the  polynomial  fit  calculations  and  the  plot. 

Y[*]  Date  and  time,  see  below. 

The  components  of  C[*]  are: 

C[l]  Number  of  files  available  on  track  6  of  the  tape. 

C[2]  Number  of  files  available  on  track  1  of  the  tape. 

CL3]  Tape  number. 

The  components  of  Y[*]  are: 

Y  (_  1  ]  Year 

Y [2]  Month 

Y[3]  Day 

Y[4]  Hour 

Y[5]  Minute 

Y[6]  Second 

Strings  employed  for  outputing  data,  in  split  precision,  to  the  tape 

are: 

A$  Contains  T[*],  H[*],  Y[*],  X[*]  and  T,  the  averaging  time. 
Special  function  keys  are  employed  as  follows: 


Key 

Action 

Purpose 

0 

*  "cont  formats" 

Restarts  programs  without  redefining 
matrices  or  loading  keys 

4 

*  sfg4 

Branches  program  so  that  other  special 
function  keys  may  be  set  and  acted  upon. 

5 

*  sfg5 

Change  order  of  polynomial. 

6 

*  sfg6 

Change  averaging  period. 

7 

*  sfg7 

Suppress  all  printing  and  plotting. 

8 

*  cfg7 

Restart  printing  and  plotting. 

9 

*  sfg9 

Print  radius  bin  center  values 

10 

*  sfglO 

Mark  a  data  tape. 
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Program  Listing: 

0:  "^>3?  \’sCoxr,L::'t-ri3  e-s  g^i'aivju'Ieg  &  ?Lor?  Ci/u"-. 

1:  iir-  A-;f  :o3l  ,r,:i3.3i '3113^01  3 

3:  :!.•  v r  i  ::i ,  uj  ,2[Jfi3],  :u*i ,7(;,i5],  ,p[o,±5i ,  :k  i 

3:  i-  ^,,1^  ,  '[.,,L3|  ,'IL'A  ,nv,±?}  ,o[i0]  ,.<f 3! 3) 

•’ :  1.13  l;3>.rfll;l ] .1*3? [2!;4>;{3|; d-.- [ 4 1 ; 42+s ( 31 ; f4*d 1 3 1 
j:  71;  »  : 

•  ■  ?*-•-»  *7  ■  ••  '/*V  11  *!  • 

-  •  ■  * —  t  '  -  *  #  • 

7:  "i;cr:jL3“:*.~.c 
i-it 

">:  foe  2  ,ei ,  l3  . 5  ,o 

I'"1:  e-t  3,  "Gate'", G ,"/" ,£2, , i!2,“  ?.t  end  of  rir." 

11:  fv.t  4,/,"' veraairvo  tire  =  ",f4.1,"  niTi-itcs" 

12:  Set  3 , "  Ftn'-o  vclteee  A  =  ",f5.3,"  volts" 

13:  cot  5,oin.2,7 

14:  t-t  7, "I sec  *-,£3,"  Pile  a", £3,"  ’’vent  i " 

15:  set  o , " Pol /no-rial  cf  order  ",f2,e 
15:  rmt  S,el5.7,s 

17:  "reset"  :cisr.  "Get  any  flans, continue. "  ;stn 
13:  if  ilalOjgto  "tancork" 

19:  ic  il35:ent  "Order  of  cclynerriaL?" , J;  iso  2 
20:  7*C 

21:  C+l-Ijcc?  All  I, G[II  ,r  (I, I| 

22:  if  tlajjor.t  "Averaging  tiie?",i':int(1.5f)*Z 
23:  dsn  "Insert  data  tare,ccntinue.”;str> 

24:  ldf  0,D^[*] 

25:  'ivc  " A inter  cn,  set  form,  continue. "  ;sco 

25:  wco  ?■’, 27,59,27 ,54 ,32,32, 32, 32,27, 77 ,27, 75,13,1 ,14 

27:  v;tr  ;  ,2/,  /3,15,tO,u,2d 

2o:  if  not  f  lef ;  2  ?-►  i4;  Lie (1.51)-*-’ 

29:  5-*<[I1  *\(2|*.*;I  ?J  ;1-»X(5| ♦aIGI  ;15-*.\l4j 

3):  1,11 ;  .|.G73*.ffl,j) ;  .3  Xa*:'  [1,3] ;  f  i  ,7 1  ;  1 1  ,f ! 

31 :  .332*;  [  l,v  I ;  .J C.>5*.;(1 ,7]  ;  .ill-*-:  [  1 j ;  .117 5-*-  [i ,  )  ;.  L>'(l,iwi 
3;:  for  .1—1 1  t-  15:7  f]  fl,.jj  ;x>t  j 

37:  :r  5--L  u  13  :  ( .  23+.  7:73 (  j  — i j )  /2*.  12,31 
3  » :  . '7 - .  2(  .-1)  ■*  ■  [3,J1  3  *-.■)•'(  j -l  )■*:  f  I :  .75*-.  77  (3-l)-»  [  5,-J ] 

75:  ii  "  13,  fl  >1; . 7:'7*"f ",  j]+.  323*',15,G| 

75:  J  ‘•(.J-l)  ’ f G i  1 ;  ,7c7*”(  >, i'lJ.733-»’('',j]  ;ro>L  J 

37:  or  T---i  to  Gofer  -;-l  to  15;  (•  [J  ,-jt-U  t-.  [I,  >1  )/o-*-  Jl,  jl  ;r.e:<t  I;  none  I 
3”:  "rr.ert":i:-tn  .3,’)-»7;  I  >:  ina  1,7,  >,£■  t  7, V ic  r  J;  :7 ;  V*  * 

’.T '  ;.fn  I  d  i.f  '=l;efa  /:)-»  ' 

,!1'ta":  ",  Id  ;4;  .•  n-Vscf:.  7 

71:  '  "  "  :ets  : ; er  —  iv.J"  :** > -1  Td  1 

-  ".res 

:7:  c!  1  *’i'(  (1| -| )  ;oli  'n'(  ):oll  '  i'(.  i;Krnf( 

-  4:  -•:-•(  ,3)  - 

.  .  J  ✓  -  •  .  -I  Is. 

d:  j  'C  l:-i  o.  2.5 ;c*  1  1 

•7:  :.r  1~-  no  '::2'C>jit  l  t:  "  ir:‘ 
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.'  C .  ,  . .  ,  j  i  •  ~  '  I  i  - 1 1  ■-  .  )  ,  •  ,  1 

1  *  :  .* ■  •*  .....  - |  ‘  ,  •  -  ii  ;  ■ :  t  i  •  i  ■ )  •  *  *  .i 

i  .  :  ll  •...  .1-^/  *u  t  1  r  ■  t  t  -•  t  '  r  •  *-  -  •  J.  f  •  '  •  ’ — 


.  .  v  /  '  L  W  •  -  >  r  *  •  1  1  *• 


IV:  V  I  .r.  :  ;t:  -j 

1  :i :  ..L-  5,  ,  ,  ir.t  (1 . 17) ,  Lit  (■'  M  ) 

IV:  it 

IV:  Vi" 

I :  i!  (cVV)  V.l: 

IV:  Mr  ;v  •,  *v  : 


17:  ic  :  '  IW  or  Ml, MM';  :fco  "  ’  m" 

V:  ];•:('  I  >. ,  l-r(  :|I,  ji)*L 

1:  x-  ,V,  M, :  •/.>,L  t(7M..),Mt(  ::/M,r  t{  :) 

:  i 1 

...  •  --  r  1  -  1 

,  i :  .»  /  r  --*•  -  * , ±  l  r 

:  r  :  1/ ; / M  , i-t{ :  V  ; f  7’ ) , ir.t (3 •  r.-. I /7) , j  it (V 

•*:,-f  ti  M:.  Mr  >1  it  :;ix  -1  to  li;  mL:  f-l>  . 

M:  3  K  yr.;x  ’ 

I  *  \\  i  .'V.  1  ’  '  *  L  j 

'  tr.  (  ! Li  •  ‘,7  *4!  l>t:  -r. 


H  i  .  ■  •  i» 

1 1  *_  •  rnll 


»  -  t  ,  f  r  i  i  ^  f  **  '  i  )  '  '  •  |  .  *  j  •  >  r-  r  '  11  \T 

..  J  V  1  i  t  I  )  *  I  .>  ..I  •  i  '  f  >  «  1  •  -t»  ■  I  t  -  i  l  l  sj  ,  l  '.V  -  •  1 


!.  ■  •:  :  ti  (/|  _I-t|  )♦'  3(  2i  l>V,3o2i-4i>) 
ii‘:  ti  (.■  f  Jf  I  |30t  ,  V4+  *2)  ;1+.W 
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150:  "rjeXT"  rrcxt  I 

151:  fts  (1) ♦A? 1 457, 460) 

152:  -l-*-{;0*Y  ;oto  "qulo" 

153:  "crv,':0-»'if;U'  X>.8;qto  ''oeldv' 

154:  for  1=1  tc  C+l;YX-KT[I]-*Y;rext  I 
155:  W  :  f  4)  ;jrro  2 

15-6:  *’bolcn”:if  a>1.  1; if  ?-:[2)>D;.*;[5]-*>Y; jr,T  2 
157 :  Mf 4)  -hi [21  * (X-.3)  *Y*\ [ 6) 

153:  if  Y<-4  cr  Y>5; jro  4 
159:  Y+6/96-Y 

150:  vtb  ;’,27,o5,int(i5«Vi)  ,int(240X)  ,int(3Y/2)  ,irit(>6i) 
161:  vtc 

152:  if  (Mr  1/1 2 i>/)  <1.5; -/to  “acv:' 

1-3:  "lulr’V.u:  ,12,13 

154 ;  of  C/.5/.1 1;  trk  0;  tvf  0 

1-3:  "flo.p":i^  .. i.:4 "rcroc" 
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2.  Plot  of  log(dN/dr)  vs.  1og(r) 


This  program  retrieves  the  data  stored  on  the  cassette  tape  by 
the  data  acquisition  program  and  produces  a  plot  of  log(dN/dr)  vs.  log(r). 
The  plot  consists  of  the  individual  data  points,  less  those  excluded,  and 
a  curve  generated  by  a  polynomial  fit  plus  end  point  extrapolations.  There 
is  no  tabulation  of  data. 

The  matrices  used  in  this  program  serve  the  same  function  as  in 
the  data  acquisition  program. 

The  following  is  a  listing  of  this  program: 


"r.V-:>2  r»X(d  /or)  FRO  5  .*u\*  Ti.Uif/': 

i:  dir  wf  ;•  ■-■! , 

-is-'  -ilfl  , ■■■I  'U 131  ,  ifl5|  ,2-r„l-j| 
i:  r;ia  l\<s)  ,.i!o,15)  ,;:,J13J 
i:  iir  rfi.i.r.1  ,  M10|  ,r(lG,10)  f  Ho] 

3:  35*',|1|  ;ill-»  ’.(21  ;4>213];0*::.I4J  ;42-tM5]  ;f,4*3ffc] 

5:  715*l;rlev 
7:  f~t  0,ir.x,r. 

1:  irt  1,0.0 
rrr.t  2,ci,£3.C,z 

Lj:  tint  3 , " Pete " ,  13 . 0 , "/" ,£2.0, "/" ,  12. 0 , 'lite",  0.0, l,:",t2.0 
11:  frt  4,/,'V'Vor>3i^  time  =  M,i'3.0,"  mirutec" 

12:  irt  6,010.2,2 

13:  f.Tt  7, "lone  File  4", 0.0,"  LVent 

14:  bnt  el vnonini  ci  order  ",l2.0,z 
15:  f^t  i),clj.7,r. 

15:  7-C 


17:  C+OIfro-  MU  ,R[I>  „•  fl,ll 


15 


23: 
^  1  . 


don  "Inrcrt  dote.  tcoc, continue '';sto 
l-’f  d,r:  ,J(*I 
ert  " ! rrcK",  i 

•iev.  "< rioter  on,  rot  Cerr, centime" ;rfc> 


vto 

'.IT 


37,5-; 

27,71 

*;Oi 


.37 ,82 ,32 ,32,32,  32 ,27 ,77 ,27 ,75 ,15 ,  J  ,14 

4  ,  •  “>  -* 

f  t  f  w.O  f  f  }  „ 

■* if  31  ;  !♦  4  f  5]  j  Ci  1 : 13+<  ( -i  1 


-’S  •',345*  fl,ll;.0375*’fl,2|;.0005*R[l,31;.094-R{l,4I;.098-*R(l,51 
25:  .102-R(l,61:.1065-Rri,71;.lll*R(l,81;.U53-Pll,9l:.12*R[l,10] 

27:  tor  J=U  to  16:Rfl  ,J-lj+.005*R[l,J|  ;rext  J 
28:  for  J=1  to  16; (.23+.025(J-l))/2*Rf2,J] 

29:  .2+.02(J-l)-*R(3,J)  ;.3+.08(J-l)-*R[4,J]  ;.25f.25(J-l)-*R[5,Jl 

30:  if  P(5,J1>1;.767*R(5,J]+.23>R[5,J1 

31:  l+(J-l)-4<[6,Jl;.767*R|6rJl+.233*Rf6,J);next  J 

32:  for  1=1  to  6;for  J=l  to  15;  (R[I,J+l|FHfI^J|  )/2*fi[I,J];next  J;next 


I 
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33:  "cycle" rent  "Begin  with  file",A;ent  “End  with  tile",Q 
34:  tor  D=A  to  B;fmt  0,10x,z;cll  'read  t ape' 

35:  for  1=1  to  & 

30:  cor  J=1  to  l5;T[I,J]-*0[  I ,J]  ;rcxc  J;next.  I 

37:  ’-'or  1=1  to  3;  cor  J=L  to-  3; (i,  M  ;ne:ct  J;r.o<t  I 

3c:  Av;|5,l]  *015,1}  ;for  J=1  tc  15;0-h.[4,J1  ;next  J 

2~i:  ,,.rat":io?  A,E;fcr  1=1  to  j;0-*-ru-*rl*r2-»r3-*-r4*r5-*-r5*r7»r'c;*r'J 

10:  tor  J=1  to  16;lV.;ir  1=6  art!  J?l;l-*r6 

41:  if  J<15;  l-n  5 

42:  if  T>1  jnc  7<5;qto  "nF-Xt" 

43:  if  1=1;-!. ■>.%;! to  "l” 

54:  if  [=6  m-l  r9 an  and  r7*t0;lag(r3./rl)*X;lo":(r7/r3)+P;l*r3;gtc  "3" 
45:  l. :>*:■: rite  "l" 

VS:  it  Of  I ,J1=0  and  r5=0:qto  "nCXt” 

47:  if  rOsOrotc  "extrao'1 


4fl:  if  0  [  I ,  J  1  =0  end  rS=l;qto  "extrao” 

40 :  Icq ( 0 1 1 , J I ) *P; loo (E  f I , J | ) -X 
50:  if  [=1  and  0<°; Jto  "3" 

51:  if  rtfqjto  "2" 

52:  a  to  ”3" 

53:  "extras"  :r  5+ l*r9:  r7+)fl  #J]+t7 ;r9F?(I,Jl*r3;  ito  "nt'Xt" 

34:  "3"  :r  -►r3;rl+.s-*rL;rxH*r2;  r 3-t-LJ<-^c3 ? r 4-»-Li*r 4 ;  if  r  5=0 : a  to 

55 :  X-IM71  ;r*t  !f|:l.  5*x 

55:  "1 " :  (r 2r2-r4rl  >  / ( rf >r2-rlrl )  -» >f  21 

57:  (rCr4-r3rl)/{rCr2-rltl)  ♦4(11;  -dll  K<*::f21-»t 

55:  -=d  to  C;:yL-^;C*-t(r,^*c  !^,P) 

53  :  P+A  ( P 1  ♦Aff.'i ;  p*X-*f-';G.6*\.. 

60:  if  (R-l  ,:%| *b  fi'-J.,r| 

61:  'l*»5;next  '< 

62:  "nfXt":next  J;rext  I 

63:  for  1=1  tc  Ji-l 

34:  for  :<=1  to  int(  (C+D/2) 

65:  if  in  and  1+602  and  I-X>0;F(  I  ,I)*F  [T-6,I+K|-Ff  (+6,1- d 
66:  u  I+;:<=>1  .md  I-.<+l>0;t'(J.,I>l)*F(l+K,I-lv+i| 

67:  if  t+r«=C  onJ  I-iOO;  Kf  I,  r +])♦?’(  £-K,I+K+l( 

63:  next  K;next  I;inv  F*F;rat  f/VC;  if  rlq7;qto  "olt" 

69:  "out":wrt  ’  ’.7" ,C[3]  ,u 

70:  wrt  "‘!.3",Y  |2|  ,Yf  3| ,Y  [1|  ,Y[4)  ,Y  |5) 

71:  wrt  ”M.4",r;wtb  10,13 

72:  cto  "no  c!n/dr" 

73:  wtfc>  M/V!Vdr";wth  M,ll),13;for  J=1  to  15;for  1=1  to  6 
74:  wrt  "*’.6",  Vll,J];mxt  I;wth  : 4,10,13;next  J;wtb  M,  10,13 
75:  "no  dn/dr": 

76:  wb  l,27,65,ir.t(-l.c5)  ,int(-120)  ,7,32,"Lcn(±S/ck)" 

77:  wto  1,27,55, l,56,7,32;wrt  ”'1.8",C;5.01-l/5*Y 
73:  tor  1=1  toC+];wth  1,27, 55, 1,56, int(3Y/2)  ,int<9GY) 

79:  wrt  "V.v"^|I|  ;Y-l/(W;next  I 

30:  wto  27, 65, irt(-3. 75)  ,int(-240)  ,0,15, "loq (radius)" 


*» 

L. 
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81:  "olt":it  f lq7;qto  "skip  print" 

82:  -1+X:  -1-»Y 

83:  wtb  M,27 ,65, int(15.</4)  ,int(24(JX)  ,0,0 

84:  if  XJO  ana  Xncdl=G;wtb  .M,"|‘',10,8,3,S;//rt  "  l.l",X;qto  +2 

v,5:  wtb 

86:  if  (X+. G5*X)  <2. 05;qto  -3 

37:  wtb  27,55,0,0, int(l. 5Y)  ,ir.t(95Y) 

83:  if  WO  and  Yrr«l=0;wrt  "M.2","-",Ynto  +2 
39:  wtb 

00:  if  (Y+.l-Y)  <5.1;atc  -3 
21:  for  1=1  to  6 
92:  for  J=1  to  15 

93:  if  r(I,J)=U  or  O|I,J]=0;qtO  "dtXV" 

94:  lai(E|I,J])*X;lcx|(0|I,Jl)-* 

95:  wtb  *',27,65,iPt(15X/4)  ,int(240X)  ,int(3Y/2)  ,int(%Y) 

96:  wto  M/’fli 

97:  "r*F.XT”:.icxt  Jjrext  I;  if  r5=0;j.r>  2 

98:  wtb  :■*, 27,55, int(15P  f7]/4)  ,int(2403(71 )  ,int(3d(81/2)  ,int(963f3| )  ,"o" 
99:  qto  "skin  crint" 

100:  ''real  tree":trk  .J;  Ldf  :),AS,5? 

131:  tor  1=1  tc  6;  for  J=1  tc  15; J+15I-16+X 
102:  stf  (-9S 1 4X+1, 4.4+4] )*r(l,J|; next  J;next  I 
103:  0-*tn;  for  1=1  tc  lo 

104 :  rtf  (AS  f  Xl+4!\,354+lX| )  -*•((  1 1  ;l+.v»4;  it  I>5;;ito  "teXi" 

105:  stf  (A$IXl+4K,384+4;<])*Y(I];l+K-*4;if  I>3;  Jto  **3X1" 

106:  itf (AS[3ol+4K,362+4Kl ) ♦X(2I-1) 

107 :  itf  (AS [ 363+4K, 364+44 1 ) +X[ 21 ] ; l+K-4 
103:  "MeXr":ne:<t  I 
105:  stf  (AS  f -157 ,460] )  ♦!’ 

110:  ret 

111:  "skip  print" :-l*X;0+Y 

112:  "crv":9-Y;U  *>.3;  Jto  "nclcri" 

113:  for  1=1  tc  C+l;YX+5fI|-Y;rRxt  I 
114:  Y- 1  {41;  T~e  3 

115:  "belch": it:  />l.l;if  •••f2]>0;‘  ffl-Y;iro  2 
116:  ‘ » f 4 1  +■* !{2 ]*(.<-. 3) ♦Y->1{6] 

117:  if  Y<-4  cr  Y>5;  lire  4 
113:  Y+6/96-Y 

119:  wtb  ::,27,;>3,int(15.'/4)  ,int(2434) ,  int(3Y/2) ,  int(%Y) 

120:  vth 

121:  if  ( X+i/1 2 0-*X ) < 1. 5 ; p to  "orv" 

122:  wto  *',12,13 
123:  next  D 
124:  uto  "c/clc" 

125:  end 
*554 
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This  program  retrieves  data  from  the  cassette  tape  for  two 
specified  files,  generates  a  polynomial  fit  to  the  data  for  each  file,  in 
log  space,  and  plots  the  difference  between  the  fits. 

The  matrices  used  in  this  program  serve  the  same  function  as  in 
the  data  acquisition  program  except  that  certain  of  them  have  been  in¬ 
creased  by  one  dimension  in  order  to  accommodate  two  data  sets.  K[*]  are 
the  coefficients  of  the  polynomial  fit  to  the  second  set  of  data. 

A  program  listing  follows: 

0:  "i To-32:  *>L?r>  U2JVS  LOGl^/Tr)  i'R:  Wm  Q:  FILE": 

1:  i.Lr.  A?  12, 4001,  i$[46f.|  ,>3[2o]  ,jf2]  ,3(0,151 ,2$  1251  ,i*[2) 

2:  .lir*  Y 1 2 ,ol  ,>  (n|  ,"[2,o,l5| , { 7 ( L(j  1 ,0(2,0,15) 

3:  Tim  l.|u)  3(6,151  ,L',0( 3j,k(10| 

4:  Jim  D[6,15]  ,46,161  ,A|101  ,F[10,10]  ,3[10|  ,1(91  ,MGj 
5:  35-2  f  1 1 ;  111-.;  f 21 ; 43-3  [3 1 ;  i-F  [4 J ;  42-:f 5 )  ;54-2  f 6 1 
G:  715-  Ijjcv  "l",.! 

7:  lTC  0  ,  lo  x ,  t- 
'3:  trnt  l,D.Q 
9:  LTit  2,cl,f5.1,z 

10:  L.rt  3, "Cote", £3.0,"/’ ,f2.3,"/' ,12.0,"  Tire", 0.0,' r2.0 ,r, 

11:  frt  4 /'Averaoirw  time  =  ",£2.0,"  and  ",£2.0,"  minutes" 

12:  fmt  6,cl0.  2,? 

13:  £;:.t  7,"Tsne  :?",£3.0,"  Files", £3.0,”  ana", £3.0,"  event  ",c25 
■14:  tint  3,"Polvrcrial  c£  order  ",f2.9,s 
15:  tmt  9,015.7,7 
15:  1*2 

17:  CH*T;rdm  />fl]  ,3fl|  ,K[I1  ,r [  1,1) 

13:  else  "Insert  Tata  tavx,ccnti.ouc";stn 

19:  ont  'Tare  -t",E(  31  ;ent  "TracK",.7 

20:  .Iso  "printer  cn,  set  tcrn,centinac" jsto 

21:  win  :»,27,6<;, 27, 04, 32,32, 32,32,27,77 ,27,70, 15, 0,14 

22:  wtb  '’,27,79 ,4 ,43,6,32 

23 :  3  »X  f  1 1  -X 1 2 )  -T [  3| ;  1-X 1 51  -.<[  f.  I;15-Xf 4 ) 

24 :  .iV?45-2(l,l| ;  .  >375-4(1 ,21 ;  .0905-: '.[  1,3 1 ;  .034-^(1,41 ;  .090-0(1,5] 

25:  . 192-' f 1 ,01 ; . 1 0i5-<(i, 71  ; .lll-n[ l,c»l ; .1155-r[l,3|  ;.12-Rfl,101 
25:  tar  .7=11  to  lo;0(l  ,J-1|  K005*:'ll,  J]  ;naxt  J 
27:  tar  J=l  to  16 ;(. 23t. 025 (0-1  ))/2 -.2(2  4) 

25:  .;+.02{J-l)-;’13,J!;.3t.G3(J-i)-f-’[4,JI;.25+.25(J-l)-rf5,J| 

23:  i i  >'15, 1] >1 ;  .7G7*k|  5,J1  +.23>r|9,J) 

3.9:  1+  { 7 -1 ) ♦: < 1 0 , J! ; . 707 | » )  + . 233 -*•  (f  5 ,  Jl ; next  J 

37  :  tor  1=1  to  <>;  lor  .7=1  to  15;  (0 1 1 ,.j+1]+G[  I,J  1  )/2-dlT,J  |;ncxt  Jjnoxt  I 

12:  "start" :ent  "i.7cnt",C?; ant  "Lacx-jrcona  file",.!  [1)  ;cnt  "Event  tile",J|21 
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33:  trk  ii;for  C=1  to  2;lu£  JfD],B$,BS;3$+AS[l'J;cll  'rort';;T:xt  D 

34 :  tor  D=1  to  2;  for  1=1  to  6 

35:  for  J=1  to  15;1*[D,I,J1-»0(l>,I  ,J];next  Jrnext  I 

36:  for  1=1  to  3;£cr  J=1  to  3;0»GfD,I,J]  ;next  J;next  I 

37:  0-*Gp,5,l]-OfD,b,l];for  J=1  co  15;0O[D,4,J) ;noxt  J;roxt  D 

33:  for  1=1  to  6; for  J=1  to  15;0(2,I,JK)[1,I  ,J]-*3[I,dl 

39:  if  SfI,J|<0;i>3[I,J| 

40:  next  J;rext  I;ror  [>1  to  2 

41:  "mat":ina  A  ,F; for  1=1  to  6;0+r0+rl-»r2*r3+r4+r5>r6*r7-*r8*r9 

42:  for  J=1  to  16;  1-^3;  if  1=6  an3  J“l;l-*-r6 

43:  if  J<15;txc  5 

44:  if  I>1  and  KSrqto  "nfAt" 

45:  if  I=l;-1.5*X;gto  "1" 

46:  if  1=6  end  r9*0  and  r7 '40 : loq ( r 8/r9 )  *X ;  log  ( r  7/r3 )  -*-P;  1+  r5 ;  g  to  "3" 
47:  1.5*X;qto  "1" 

43:  if  O(C,I,J]=0  and  r6=0;gto  "n£Xt" 

49:  if  r9tf0;qto  "extrao" 

50:  if  0[D,I,J]=0  and  r6=l;qto  "extrao" 

51:  lcxi(0|DfI,Jl)-*P;loci(E[I,Jl)-X 
52:  if  1=1  and  J<9;qto  "3" 

53:  if  I#6;qto  "2" 

54:  qto  "3" 

55:  "extrao" :r9fl-r9; r7+on,I,J]-r7;r8+efI ,JW8;gto  "next" 

56:  "3"  :r0+XX+t0;  rl+/.-*rl; r2+l-»r2; r3+PX*r3;r4+P*t4; if  r5=0;ato  "2" 

57:  X-nH71  ;p*Bf31:1.5*X 

53:  ”1 " :  (r 3 r2-r1 rl > / ( r0r2-rlrl )  *>1(2} 

53:  (rur4-r3rl)/(r0r2-rlrl)  •*'•'(  11  ;'i fl) +/*'’.[ 2) 

60:  "2":for  X=0  to  C;C+l-X*3;Gtt’tF.,<] *F (!?,p| 

61:  P+A  (R)  >A[r>) ;  P*>-*p;  ZX~: 

•52:  if  KrCjC+FfO-l  ,!'|  *F|  t'-l  ,R| 

53:  ;  next  X 

64:  "n£Xt"  :ac>:t  J;r.oxt  I 

05:  ter  1=1  to  C+l 

55:  for  X=i  to  ir.t{  (C*-l)/2) 

*57:  it  141  and  I+KC3+2  anil  I-K>0;  f  [  1 ,1 1  ♦f  f  I-X,  I-H1]  »F [  I+K,  I-K] 

53:  if  I+K«tl  rr>;i  I-;;+l>0;FlI,I+H*F( ItK,I-X+l] 

63 :  if  I+K<C  anl  I-X>0;  F  [  I ,  I+-1J  -*F[  l-X ,  I+X+l] 

70:  next  K;next  I;inv  F+F;i£  Ol;rat  PAC;j.:r>  2 

71:  rat  FVK 

72:  next  C;aro  KH>3 

73:  "out" :or t  "1\7",C[31  ,J|21  ,7[1]  ,CS 

74 :  Vf 2 ,4 1  — Y [1,4]  +2 jif f 2,5 J-Y f l,51+5d*i2*Z 

75:  wtt  "■•1.3,’,Yf2,2i,Yf2,3|  ,Y{2,1]  ,*[2,41  ,Y  [2,5] 

74:  fat  5,"  dT=",f3.0,"  .ninutes" rwrt 
7?:  wrt  "”.4",5[21,3[ll;vt^  1,10,13 
73:  ntc  "no  in/dr" 

79:  wtb  'V'drVdr" ;wtb  1,30, 13;  for  J=1  to  15;for  1=1  to  6 
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30:  wrt  ‘'M.5",T[I ,J] ;next  I;wtb  '4,10,13;next  J;wtb  M,1Q,13 
81:  "ro  dn/dr": 

82:  wtb  M  ,27 , 55  ,-4 ,0 , 7 ,32 , "del ta [ log  ( dfl/dr ) ) " 

33:  wtb  M,27 ,65,1,56, 7,32;wrt  ,,M.8",C;5.01-1/6*Y 
34:  for  1=1  to  C+l;wtb  M,27,55,l,56,int(3i/2)  ,int(96Y) 

85:  wrt  "M.9” ,GfI) ;Y-l/5-Y;nGxt  I 

36:  wtb  M,27,65,~4 ,-210 ,-1,-16, "log(radius) " 

87:  "olt":ic  flq7;qto  "skio  nrint" 

88:  -1-*X;-4*Y 

89:  wtb  ?!,27,65,int(15,V4)  ,int(24QX)  ,0,0 

90:  if  X40  ana  <rcr!l=0;ivto  ;',"i",10,3,3,3;,'/rt  ’M.l",X;nto  +2 

91:  wtb  n,"-" 

92:  if  (X+.05*X) <2.05;gto  -3 

93:  wtb  27,55, 0,C, int (1. 5Y)  ,int(S5Y) 

54:  if  YtO  anc  Yrcrl=0;wrt  "il.2","-",.5*Y;qto  +2 
95:  wtb  N,-r 

96:  if  (Y+.l-Y)  <5.1;gto  -3 
q7:  gto  "nkio  rrirt" 

°8:  "rcrt" : for  1=1  to  6; for  3=1  to  15;J+15I-16+i\ 

99:  str  (A8I’3,4<t-l,4:<+41)  -»I  (0,1  ,J]  :naxt  J;ncxt  I 
100:  0*X;for  1=1  to  16 

101:  stf  (AS  f.5, 361+4X , 364  *-4K) )  -*i i  f I  ] ;  1+K+X;  if  I>5;qtO  ’McXT" 

102:  st£  (A$(1',361+4n,364+4K]  )-»Y  [0,1]  ;1+X-*K;  if  l>3;j|tc  "IdoXT" 

103:  itf  ('1$  [0, 251+4K , 352+4:4] ) ♦X (21-1 1 

104:  itf  (A$|D,:r>>4K,364+4K]  )-Xf2I]  jl+K-X 

105:  "“aXrrroxc  I 

103:  Gtf  (?,$  (0,457 ,460] )  -2  (D) 

107:  rot 

103:  "bKio  or  int" :-l*X;  0+Y 
100:  ”crv" 

110:  fcr  1=1  co  C+l; YX-Ki(I]-»Y;re:<t  I 
111:  if  Y<-4  or  Y>5;i^o  4 
112:  ?*(Y +6/03)  v/ 

113:  wtb  :!,27,j5,int(15V4),int(24(K),int(3Y/2)  ,int(%Y) 

114:  wtb  m,*.« 

115:  "nolc"  :if  CXH/12>K)<.35nto  "orv" 

116:  wtb  ‘'1,12,13 
117:  -jtc  "stsrt" 

118:  end 
*30256 
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4.  Tabulation  of  A(1og(dN/dr) )  vs.  r: 

This  program  retrieves  data  from  the  cassette  tape  for  two 
specified  files  and  tabulates  differences  between  them  in  two  formats.  In 
the  first,  a  bin  by  bin  difference  between  log(dN/dr)  for  the  two  files  is 
tabulated  vs.  r.  In  the  second,  a  tabulation  of  the  differences  between 
polynomial  fits  (as  described  in  the  preceding  program)  are  given.  The 
matrices  serve  the  same  functions  as  described  above. 

A  program  listing  follows: 

C:  ".i  XT-32:  IVJVTSS  VUA  (tOZ(S'ViX) )  v3.  r  VC  T::~  v\WC3  CP  TdE": 

I:  "HTLW  r.'r,  713  V;  iriin  FVC 1X3 ;  iVi  VAPI^  RAC  II" : 

2:  ■:!-  A5(2,  ;*.0)  ,2Sf4£.f]  ,tSf3S|  ,/f2|  ,rfC,15)  ,?$[7S)  ,112] 

>.  ..L".  V[2,5I  ,;:fG|,if:,o,i5I,  :fl6|,:{2,5,151 
"12)  /•.rr.,13.1  31 ,  nTICI  ,Lv(23! 

3:  :i-  nfC,i-31  ,>[6,131  ,M10l ,2110,101  ,3(101  ,I(9|  ,'(6| 


701*1;. lev  "  ", 
.at  C, lux,: 


> :  : 

tot  1,13.0 

15: 

fmt 

2,  cl,  14. 

.5,2 

11: 

3,“bV 

fO.df  / 

cy  r  n  /ii  r\ 

;  i.4ii j  f  /  r  Lz~.  U 

,"  Tiro" ,  £3 . 0 , " : " ,  f2.  J , 2 

12: 

bt 

6,011.2 

9  /- 

13: 

rot 

7, '"iano 

4ii  r\  u 

t  f Lj, J  r 

Tiler", £3.0," 

and", f3.0,"  I.Vcnt  ",c25 

14: 

fmt 

9,ol5.7 

15: 

7*3 

15: 

C+l-> 

•  I ;  r-h  A 

(1 1  »o(I) 

17: 

dso 

"Insert 

data  tare 

,ccntinuc";stD 

18: 

ent 

'Taco  tl 

",C(3| ;ent 

"Track”,  J 

19: 

roc 

"Printer  on,  set 

form, continue" 

;sto 

20: 

wtb 

M,27,69 

,27,34,32, 

32,32,32,27,77 

,27,76,15,0  ,14 

21: 

wtb 

T, 27,79 

,4,43,5,32 

22 :  . 0345*R { 1 , 1] ;  .087  3*R  [  1 , 2] ; . 0905*» [1,31;.  CW4*Rf  1,41 ;  .093-R [1,51 
23:  .  102*61  1,61 ;  .1CS5*R[1,71 ;  .111*0(1, >31 ;  .1155*R(1,2]  ;.12*R[1,10[ 

24 r  for  J=.ll  tc  lS;P.fl,J-lJ+.005*S[l,JJ;next  J 

25:  tor  7=1  to  16 ;  ( .  22+.025 (J-l) )  /?+'.>. [2  ,Jl 

26:  .  2-f  .02  (J-l )  <■’.>]  3,  J[ ;  .2t.08(  J-l)  *0  ( 4 ,  J 1 ;  .25+ .  23(1-1)  *r  f  3,  J1 

27:  it  r{3,J?>l;.?67*7[5,Jl+.233*~[5,JJ 

25:  1+(J-1)*  '15, J( ; .  7 57*9  [  5 ,  T[ + .  233*  ■’  f  5 ,  J 1  ,-noxt  J 

2°:  icr  1=1  to  5:  for  1=1  to  15;  (^  f  T ,  7-+T  J  +-rr(  r ,  1}  )/2V’fI,  Iffiioxt  J;noxt  I 

I):  /=!  to  2;..tc  ,10,10 

31:  one  ”  „ mV, 35;  25* "loll  ,2i-lon(CX)|  *r> 

r:  -r,t  •*  w -rr  "wu  hvvip.i 
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trk  afcr  >1  to  2;l’f  jr;]  ;cll  'sort'; next 

for  r-1  co  2:  for  I  -l  to  5 

for  J=1  cc  lb;  rfi:-,l , J 1  *0(0,1 ,J]  ;rcxt  J ;ro:<c  I 

cor  I-L  co  3; fee  .1=1  to  3;.>2>[f  ,I,Jl  ;r.cxt  J;.xxt  I 

l.o ,5 ,11  *.:  i ‘-.G  ,11 ; for  J=1  to  lb;!>OfLy:,J|;ne.<c  J;noxt  C 
tor  1=1  co  5; for  J=1  to  13;if  J[2,I,J]=0  cr  0[3,I,J)=0;-1000*3[ t,J|;inr>  2 
iaiPf  2, 1  rj  I )  -logo  II  ,1 ,:  I  ? +:[!,. j] 
next  Jjnoxt  I; tor  D=1  tc  2 

"nat":ina  A , F ; for  1=1  tc  i5;O*r0+rl*r2*r3*r4*r5*r5*r7*r3*r& 

for  .J=l  to  15;  1*3;  if  1=5  and  Jifl;l*rG 

if  J<15;  t.'.r  3 

if  I>1  end  I <6; ate  "rEAt" 

if  1=1  2 ;gtc  "1" 

it  1=1  mo  r.l?0  anu  r7  „l  0 ; Ioq ( r 2. /r> )  *\ ; Ion ( r 7/r  1 )  *P;  1* rj ; ^ tc  "3" 

it  jfd,I,Jl=.)  art’.  r3=o;  3co  "rJIiit'' 

it  r')iC;.;co  “cxtrtr" 

it  ■  if  j,I  ,j  1=5  art:  ca=l;  jtc  "over j 

lc«.t(  r '  ,  T  ,  jl  )♦  !v;lcx.(f  |1,J|  )  :x 

il  1=1  an..  JC>;  ttc  '’3" 

if  Ifvjvjto 

ttC  "J" 

“c\'tra’'.-’;r.H-i*rc; r 7 K ; Id:,/,  '*t7;rdH-[I,J.j*r3;?tc  "n'iAt" 

'*3"  trO Kv>rj;r]K*rl  •  r.'''  t*t ::2f  A*r3;r4+f*r4;if  r5=0 ;gto  "2" 

;>2  17)  ;>•  ft] 

"1“ :  (r2ri-r4 rl )  /( r Jr2-rlri)  *  :[  2  i 
( rOr-t-  rj  rl )  /( r.'r2-ri  rl )  -  if  1 1 ;  :  II !  +.*.*:!  21  ♦+> 

•'2"tfor  tc  ) 

•+.'(o^-.f.'];2V*?;.:<*'3 

if  <"C;:.4v'f;-l/'l*«M^i#r;1 

T-*S'*5;  nay  t  1 

t  r 

tar  [=i  to  2+1 
far  >1  f:  i.ot(  (2+D/21 

it  T  ■  l  and  r+t'O+2  an.1  J,ll*i  [l-2,I+a'l*Hl+  :,I-'  l 

if  rv’«+l  yJ  I-..+ IbCsP n, I+l|*7! T+'^l-U+H 

it  I  +••'<— 3  art:  I- 2>Cit-’|I *ff  T- C,  1+2+1) 
neat  :.;na:-.t  Ijinv  F*d;il  >i;..oc  <:7*0? iro  I 
re  t  i- 

02/t  3;or?  .\-G*3 

"oct":./rL  ">;.?",C[31  ,J[2I  ,J[1|  ,25 

'/ft,  I |  -'/ ( 1  /•  1  *Z ;x  f 2 , 5 1 -r  r  1 , 3 1  +3C*3*Z 

•art  '-..2" ,*[2,2]  ,r(2,3|  ,Y[2,1]  ,*[2,4] ,/ [2,5] 

Cr t  5,"  51="  ,0.0,"  ..rinutet" ;vrt  ";.5",5 
•/to  i  ,10,13 

v/to  I • , "f ol to ( log ./dr ) ) " ; *.;to  ",]0,13;£or  J=2  tc  15; for  1=1  tc  j 
it:  i=4;  :co  "  lT,.Ay" 
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SO:  iu  3(1,. Jk-J?‘3;wtn  i,"  ",42,"  ";j.rc  2 

SI:  wrt  "”.o",3U,J| 

32:  "CV'L >'V"  tncxt  Xjwtt*  10,13: next  Jjvtb  3,10,13 
S3:  -no  'A.'V* 

34:  ''rcrt":Lcr  1-1  to  j:  ter  J=1  tc  l'3;0+13l-r5-*-< 

-jS:  -tc(.!'Sf3,4.<+lf-».<+4))-*0[O,I,Jl:roxt  J;next  1 
J5:  t'+.’rfer  1-3  trio 

17:  st£  3S1K,<, 354+4*;] 3  *i|I)  ;3  +  .;*K;if  I»;cxo  "STcxr" 

S3:  r.a  ('k.(3, 331+44, ZSlf  kn+sT-', II  I>3:?to  "  IqAv" 

30:  itf  (•■'■:  (0 ,  Jo  1+4.4 , 352+4X  ] )  -X  f 21-1] 

00:  itf  (AS  |0,363+4 <,354  +4*l)-»>* {211  ;l+K-*-X 

91:  "  JeXT"noxt  I 

32 :  r.tf  (AS ( 0  A 57 ,460] )  *2- f 01 

03:  rot 

94:  "'JPV*  :for  D=1  tc 'l;£cr  J~1  to  5;  .05(S(  J-1)+3)-1.05*a;0*^ 

95:  for  1=1  tc  3+l;tXKSfI !  *f;  next  I 

96:  i'*3f J,3 1 

97:  next  J;noxt  0 

93:  v;tb  ”,1 3, "Fitted  CtTV?";wtb  44.10 ,13:f car  J=1  to3:for  1=1  to  5 

99:  if  3[I,Jl=0;wtb  M,"  ,,,42,"  " ; "no  2 

100:  '•Art  H’1.6",sn,Jl 

101:  next  I;wtf>  4, 10, 13; next  J 

102:  next  V 

103:  wtb  f\13 

104:  a  to  " SlAFl " 

105:  end 
*6233 
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5.  Plot  log  (A(dN/dr))  vs.  loq(r) 


This  program  retrieves  data  from  two  specified  files  on  the  cas¬ 
sette  data  tape  and  generates  a  polynomial  fit,  in  log  space,  for  each  of 
the  two  files.  It  then  plots,  point  by  point,  the  logarithm  of  the 
difference  between  the  antilogs  of  these  two  functions  versus  log(r).  The 
matrices  serve  the  same  functions  as  described  above. 

A  program  listing  follows: 


0: 

1: 

"V  . 

3: 

4: 

5: 

5: 

7: 

6: 

10: 
11: 
12: 
13: 
14: 
15: 
16 : 


";>\3-32:  PfOI'J  LOG(3£L'I7'  (dH/dr) )  FPCM  TATA  ZW  FIlC": 

din  ,■'•$[2,4621 ,33[4GC)  ,u$[3f>]  ,1(2]  ,$[6,151  ,CSf25]  ,$[21 

d im  VI2 ,01  ,3  1  fil , k’ f 2 ,5 ,15)  /illol  ,01 2,6,15) 

dim  481  ,$[6,151  ,$,$[31 ,3 [101 

'in  D ( 6 , 15 1  ,2(6,161  ,\flO)  ,P|  10,101  ,$[10]  ,IfS|  ,  -1(6j 

25*?.  n  1 ;  1 ll-:)  [  2 1 ;  4 5-3  ( 3 1 ;  OP  (4  ] ;  42-2 f  5)  ;64-S (61 

715-*v’;3cv 

fmt  0,10.x, z 

fmt  1 ,  f3. 3 

tnt  2, cl, t 4. 0,z 


t'nt  3,"Dote",f3.0,"/,‘,  12.0,"/',  12.0,"  Tine",£3.G,":",f2.0 
rnt  4 ,'Vveraqinq  tine  =  ",12.0,"  and  '',£2.0,"  minutes" 


tnt  6,el0.2,z 

rat  7,"Tnoe  if", 0.0,"  Files", 0.0,"  and", 0.0,"  Event  ",c25 
rmt  a/'Pclynaaial  of  order  ",f2.0,z 


imt  9,el3.7,z 


7-C 


17:  3+l-I;rdm  *MI1  ,G|I]  ,K[Il  ,F[I,I] 

13:  dsn  "Insert  data  tace,ccntinus";sto 
19:  snt  "Taco  .?",Cf  31  ;cnt  "Track",./ 

20:  dsn  "rt inter  on,  set  form, continue" :sto 

21:  wtb  i’,27, 59,27,34,32, 32, 32, 32, 27, 77 ,27, 76,13,0 ,14 

22:  wtb  *:,27,79 ,4 ,48,6,32 

23 :  >X  f  1 1  -a [  21  -<  f  31 ;  1* If  51  [ 5  ] ;  15-'<  (4 1 

24 :  . 084 3-2 (1,11;. 037  3*0 (1,2 1 ; .0305-5(1,31 ; .034-3(1,41 ; . 096-3 (1,51 

25 :  .  102 *? ( 1,6 1 : .  1 055-3  [  1 , 7]  ;  .111-3  f  1 ,3] ; .  1155-5 (1,31;.  12-P  [  1, 10 J 

25:  for  0=11  to  15;';(1  ,.7-ll  f-. 605-3(1, 31  ;ne.\t  J 

27:  for  )=1  co  16; (.:3f-.025(J-l))/2-x[2rJl 

23:  .f+.<>2(  1-1  )-3[3,Jl ;  ,3f-.03(  f-l)-3[4  ,J] ;  .25'-.25(  J-l)*?(  5,3) 

29:  if  i' (5 ,  J 1  >1 ; .  757*P ( 5 ,  J 1  +-.  233-5  [5 ,  J1 

30:  1+(J-1)-P[6,J|  ;.7>7*::(5rj  |+.233-rtfb,J]  ;next  J 

31:  for  t=l  to  5;ior  3=1  tc  15;  (»’.[!, Jt-ll+Kf I  ,i]  )/2-i$(I,Jl  ;next  J;next  I 
32:  "rtert":ont  "iver.t",$S;ent  "lackarcurvl  file" ,j  |lj ;ont  "Event  £ilc",.J(2] 
33:  trk  ><;for  3=1  co  2;luf  J  [i)|  ,53,bS;3$-/$(o)  ;cll  'sort';next  L> 
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34:  for  0=1  to  2; for  1=1  to  5 

35:  for  3=1  to  15;T[D,I  ,Jj-*O[0,I,J]  ;next  J;next  I 

36:  for  1=1  to  3;for  J=1  to  3;Q*0[D,I  ,J]  ;next  J;rext  I 

37:  0-0(0 ,5,11*0(0, 6, l];for  J=1  to  15;0*O[D,4,J];next  Jrnext  D 

38:  for  1=1  to6;for  J=1  to  15;0{2,I,J]-0[1,I,J]-*3(I,J) 

39:  if  S[I,J ]<0;0*S[I,J1 
40:  next  J;next  I;for  D=1  to  2 

41:  "TOt"  :ina  A,F;for  1=1  to  6;G*r0*rl*r2*r3*r4*t5*r6-*t7-*-r3*r3 

42:  for  J=1  to  IS; 1*3; if  1=6  and  -J*l;l*r6 

43:  if  J<15;  jrr  5 

44:  if  I>1  and  KG;  a  to  "rfiXt" 

45:  if  I=l;-1.5*X;qto  "l" 

46:  if  1=6  and  r9s0  and  r7  •*  0 ; Ion ( r 8/r9 ) *X ; log (r 7/r 3 )  *P;  1*  r5 ;  g to  "3" 
47:  1.5-X;gto  "1“ 

48:  if  G[c,I,3)=0  and  r‘6=0;gto  "nRXt" 

49:  it  r5W;gto  "cxtran" 

50:  if  Q[D,I,J1=0  arrl  r6=l;qto  "extrao" 

51 :  Icq Q (0, I ,J|)* F;loa Of  (1 , 3] ) *.< 

52:  if  1=1  and  JO;qtc  "3" 

33:  if  K6;cto  "2" 

54:  qto  "2" 

55:  "extrao":r3+l*r9;r7tO[9,I,j]*r7 ;rb-Hf  (I,  Jl*r&;4to  "nCXt" 

56:  "3"  :rO+>0(*rQ;tHX*ri;r2+i*r2;r3+PX*r3;r4+]f*r4;if  r5=0;3to  "2” 

57:  X*P(7] ;P*3f81  ;1.5*X 

58:  "1": (r3r2-r4rl)/(r0r2-rlrl }*M(2J 

53:  ( r Or 4- r3 rl l / ( rQ r2-rl rl )  ♦'  1  [  1] ;'1fl]+A*M[21*P 

60:  "2": for  X=0  to  C;C+l-X-^;G+i?U'»P]-FlP,Fl 

61:  P  t-A[P.l  *A  [F] ;  P*X*P;GX*C 

62:  if  uip2;'3+F[P-l,Pl-F[F-l,P) 

63:  G*X*G;ncxt  .< 

64:  "n2Xt“:ncxt  J;naxt  I 

65:  for  1=1  to  Ctl 

55:  for  <=1  to  int((Ctl)/2) 

57:  if  Irl  and  I+X02  and  I-K>9;  F[  I,Il*f(I-X,  £4K)*P[I+X,I-r:l 
63:  if  I+K«+l  end  I-X+l>d;F(I,I+l}*Fft+K,I-X+l) 

59:  if  I+K«  and  I-^>0;F[I,I+1]*F(I-X,I+K+1] 

70:  next  X ; next  I;  inv  F*F;  if  D=l;oat  FA*J  ;i,To  2 
71:  rot  £7\*:< 

72:  next  0 

73:  "oafr./rt  "X.7"  ,2(31  ,J[2]  ,J  [1]  ,3$ 

74:  f [2,41 -i’ll, •*] *2;y[2,5)-fl  1,51+60*2-3 

75:  vrt  "r  .3" ,f  [2,21  ,Y[2,31  ,Y[2,U  ,Y(2,-I|  ,/{2,5] 

76:  frnt  5,"  J'l=",f3.0,"  minutes" ;vrt  "’1.5",X 
77:  art  •"•■.4',,742|,3fl);.vtb  1,19,13 
7P:  oto  "no  dn/.ir” 

79:  w to  '’,'V’:;/dr"(n/tt)  1,10,13; for  J=1  to  15; for  1=1  to  6 
•39:  v.Tt  ",,.5",T f I ,J1  ;next  I ; -.vtto  ’1,10, 13;. next  J;wtb  M,10,13 
31;  "no  dn/dr"; 
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32:  wtb  v1, 27, 65  ,-4,0,7,32 ,"loa [delta(df/dr)  1 " 

33:  wtb  M,27,65 ,-4,-210 ,-1,-15,  "log(radius) " 

84:  "plt":if  flq7;qtc  "skio  print" 

85:  -1*X;-4*Y 

86:  wtb  y, 27,65, int(15X/4) ,int(240X) ,0,0 

37:  if  XSO  and  Xnafl=0;wtb  rVT',10,3,3,3;wrt  “M.l",X;gto  +2 
33:  web  .* 

39:  if  (X+.  05*X)  <2 . C5;  a  to  -3 

00:  wtb  ■',27, 65,0,0, irt(l. 5Y)  ,int(96Y) 

SI:  if  Y?0  and  Yrrocl=0;v/rt  ,^.2",,'-",Y;jto  +2 
92:  wtb  r,,,r 

93:  if  (Y+.l-Y)  <5.1; a to  -3 
94:  qto  "skip  print" 

95:  "rcrt'bfor  1=1  to  5;fcr  J=1  tc  15 ;Jtl5I-16-*K 
96:  ntf  (AS[0,4K+l,4;:+4))*l-p,I,J]  ;next  J;next  I 
97:  tor  1=1  to  Lo 

98:  stf(^[0,361+4X,3o4+4K|)-*<f(I|;l+.<-H<;if  I>5;.jto  eXi" 

39:  ntf  (AS [0,251++<,354+4Xl )  ♦Y[0,J)  ;i+!>-*-.\;i£  I>3;jto  "  »cXr" 

100:  i  t£  ( [0, 351  M.x ,  362  w;.\ ))■*■.<  f  2 1  -1 1 

101 :  itf  (A3  (C,  353+4K , 354-KK] )  -X  [21 1 ;  1+K-X 

102:  "NeXT" : next  I 

103:  ctf (AS [0,457,450] )*ZfCj 

104:  ret 

105:  "skio  orint":-l*X 
106:  "crv":0-*i-»Z 

107:  for  1=1  to  C+l;YX+G[Il-^;2X+i<[I]-2;next  I 
103:  tn"2-tn*i*2;if  Z<=0;Jto  "nelo" 

109:  loa(Z) *i 

110:  if  Y<-4  or  Y>5;iniD  4 

111:  Y+6/96-Y 

112:  wtb  M,27,55,int(15X/4) ,int(240X) ,int(3Y/2)  ,int(96Y) 
113:  wtb 

114:  "helo" :if  (X+l/12C*<)<.3Mte  "erv" 

115:  wtb  ",12,13 
116:  qto  "start" 

117:  r?rrl 
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6.  Tabulation  of  log  (A(dN/dr))  vs.  r 


This  program  takes  data  from  two  specified  files  on  the  cassette 
tape  and  tabulates  differences  between  them  in  two  formats.  In  the  first, 
the  logarithm  of  the  bin  by  bin  differences  between  dN/dr  is  tabulated.  In 
the  second,  a  polynomial  fit  of  log(dN/dr)  vs.  log(r)  is  generated  for 
each  of  the  two  data  sets  and  the  logarithm  of  the  differences  of  dN/dr, 
calculated  from  these  two  functions  at  specified  radius  values,  is 
tabulated.  The  matrices  employed  are  the  same  as  in  the  previous  programs. 

A  program  listing  follows: 


f»:  LXCGCAC/cr) )  vf>.  r  mi  VMXis:5  u?" : 

I:  "2'E  HOC*  T.-  T’I.  FTi*  10  7. US  flCiTG:  Ljj'  W.1ZU5  PAH  I" : 

?:  Cir.’  ’$[?.,*- 1,  ,>*[361 ,1  [21  ,£[<>,35)  ,27(231  rZ\2] 

3:  ai~.  x\2, 61  ,*M,ll2, 5,151  ,  L61 , 3(2,3,151 
h:  Ji-  ~foi 

5:  .\L-  5(6,151  ,  (..,131  ( 1,-1  ,'2!  13,10)  ,J{1C]  JC1  ,  TGI 

o:  “  "*i.‘5 

7:  35*0  111  ;1  LI*.- Cl ;  43*- f 3 1 ; >  *  a  [ 4 1 :  ;2* ;.  1 5 1 : 54*... ( o ] 

5:  701*1  ;.icv  *' 

0:  L~t  G,lCx,r 


10: 

Let 

J. ,  L  J  •  U 

11: 

Lot 

2  ,cl ,  •  0  fz 

12: 

tmt 

3,"::atc",f3.0 

,  L2 . 0 , 

r2.0, "  Tire" ,  i'3. 0 , ": ", t2 

.o,Z 

13: 

rent 

o  ,ell.  1 

Li: 

t'rr.t 

7, "[are  ?",f3. 

•  0," 

Files 

", f3.G,"  eixi", f3.C,"  [Vent  ” 

,c25 

15: 

irt 

iV'tolyncrial 

ot 

oroer 

", £2.0,z 

16: 

rnt 

9,el5.7 ,z 

17: 

7*2 

IS: 

10: 

20: 

21: 


C+l*I;rdm  A(I1  ,3[l]  ,K[H  ,P[I,Il 
iso  "Insert  Oata  tacc,ccntirua";st) 
ont  "Tree  s",2[31;e-.t  "Track",  * 

>'?o  'Tr  inter  cn,  set  f erm, continue" ;cto 

wtb  .'>27, 10 ,27,34,32 ,32, 32,32,27, 77 ,27, 75,15,1), 14 

wtr  ", 27,??, 4, 43 ,»i,32 

>.< f  1 1  *< ( 21  *•:( 2] ;  1*< 1 51  *'.[ 0 1 ; l'J*<[4 1 

.  J!  !-13*e  f  l ,  i  | ; .  00 7  3*  ’  f  1 ,2 ) ; .  0  .  <}3*  '  ( 1, 3 1 ; .  30  4*>T  1 , 1 1  ; .  >  [  1 , 51 

.  102**  1 1  ,.51 ; .  100.5*  ’  f  L  ,7 ! :  .111*  U  L  ,n] ; .  11 35*T  1 ,01  ,* .  12*'  11,  Cl 

Lvr  T— 13  to  13;Pri,-.)~ll  KiJ05*7fl,jJ;rt!xt  J 

for  J=1  VO  14;  325(1-1  ))/2*.’f2,Jl 

.2-*-.,.? ( J~1  )**.( 3,  ;1  ;.3+.0L(.J-L)*i  f-i  ,J  1 ;  .25+. 23(0-1  )*,:[  3,J) 

if  .'f5,  j  I>1;  .7.47*'  1  5,JH.223*.’!5,-I| 

1*  (i-1 )  *  ’  ( 5  ,J| ; . 7T7*f :  I  o ,  J I  V  233*4  f  5  ,J  1  jraxt  J 


FHbCU»l*C  bU4»K-NOr  FlU«JbD 
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32:  for  1=1  to  5;tor  J=1  to  13;  (P  ( I ,  J+l]  +r:[i  ,J  ]  )/2-*3[l,J];next  J;ne>ct  I 

33:  "OIARi"'  :i'or  \/=l  tc  2;'., to  >’,10,10 

34:  ont  "FVcnt",C3;CSsJ$|l,24-len(C5)l <5 

35:  ent  "iirtCKSPJUi-D  FID:",  J[lj  ;ont  "tVENT  PILE",  J 12] 

36:  trk  Wjt'or  '>1  to  2;ldf  J[D],G$,E$;B$-*/'$[D] ;cll  'sort';rext  D 

37:  for  0=1  to  2;  tor  1=1  to  6 

33:  for  J=1  to  15;  T(D,I,  J|-*G{D,I,J]  ;ncxt  J;next  I 

39:  for  1=1  to  3; for  J=1  to  3;0*O[C,I,J) ;next  J;rext  I 

40:  0-*C[D,5,l]-*G(D,6ll]  ;for  J=1  to  15;0*OfP,4  ,J)  ;next  J;next  D 

41:  for  1=1  to  6; tor  J=1  to  13;0[2, I,Jj-0[l,I,J]-*S fI,J] 

42:  if  SfI,J)<=C;C-*3[I,J];jrro  2 
43:  lcq(SlI,J])*3|I,J) 

44:  next  J;ncxt  I;for  !>1  to  2 

45:  "rat"  :im  A , p; for  1=1  tc  5  ;0-*rr)-*rl*r2-*r  3*r4 -*r5*r 3*  r7 -* «-* r9 
46:  for  3=1  to  lri;l-»j;if  i=5  aal  J‘i;l*rb 
47:  if  J<16;yin  5 


13:  if  I>1  and  i<6;^to 
4D:  if  I=l;-L.>.4;otc  "l" 

50:  if  I=fc  cnc.  r)$6  cru  r7 i J; Lo.? ( ro/r )) -*rf;lc:j (c 7/c3 ) +i~, l-r5 ; 9 to  "3" 
51:  l.5-*-.v,v<to  "1" 

32:  it  •’>(0,1  aai  r6-:);  ;to  "a’/t" 

53:  if  rJ-tJ;  to  "caLtoo" 

34:  if  3(3,1 ,0 !-  3  ort>  rO=i;tto  "oxtrno" 

55:  loq(f  f.3,I,J!  )*;  ;loi(t41  ,  >)  )-*■! 

3o:  if  1=1  ana  JOa.to  "3"  v 

57:  if  1,0; -to  "2" 

33:  e  to  "5" 


>4:  "oxtr<vV*  ;r  •  •+  l*r  4;  c?+7(  <,I,Jl-*r7;roH  ( (,.Jj*ro;3to  "nf.xt" 
ml:  "3"  : rt) k r< r2<-l-»r2;r3t.‘;x>r3;r4*-H*r4;i!!  r5=2;.jto 
61 :  y-*-j  [7] ;  ( 6 1 ;  1 . 5*3 

52:  "1": (r3r2-r 4rl)/(rOr2-rlrl) -*1(2| 

63 :  ( r0r4-r3 rl ) /( r0r2-rl  rl)  v-;(  1 ) ;  •:  1 1 )  +x*:’(  21  ♦  P 
64:  "2": for  <=0  to  C; C OF 1P,F| [ R , P) 

65:  PM(R1M(F|;P*X*P;GX-C 

56:  if  K<fC;G+r’(R~l  ,R)*P(P-1,HJ 

67:  G*X*G;next  K 

68:  "nEXt":noxt  J;next  I 

69:  for  1=1  to  Ol 

70:  for  K=1  to  int((C+l)/2) 

71:  if  lid  and  I+K<C+2  ana  I-K>0;f  l  I,Il*t[I-K,I+Kl-*F(I+K,I-K] 
72:  if  I+K<=C+1  cni  I-K+1>U;P( 1 ,1+11  *F[ I+K,I-K+1| 

73:  if  IH«C  and  I-K>0;Pf  I,I+1)-*F( I-K,I+K+1J 
74:  next  K;next  I;inv  F+F; if  D=l;rrat  FA*G;-jiTn  2 
75:  mat  FA*K 
76:  next  D 

77:  "out"n;rt  "V.7",C(3!  ,J{2|  ,J[ll  rC$ 

78:  Yf2,4|-Y(l,4|*Z;f(2,51-Y( 1,5|+60*Z+Z 

79:  wrt  "M.3\YU,2|  ,Y(2,31  ,Y[2,li  ,Y12,4|  ,Y[2,5} 


11  All 

4. 
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90:  fmt  5,"  rfT=",f3.0,"  minutes" ;wrt  "M.5",Z 
81:  wtto  M,10,13 

82:  wtb  Mf"lrn|delta(JN/dr) ) ";wtio  M,10,13;for  J=2  to  15;for  1=1  to  6 
S3:  if  I=4;ato  "CNWAY" 

84:  if  S[I,.T)=0;wtb  M,"  ",42,"  ";imo  2 

85:  wrt  ,,:\o",3H,Jl 

35:  "GAL  AY"  :r»cxt  T;wto  1,10, 13; next  J;wtb  .'•■,10,13 
37:  nto  "2PV" 

OR:  "cert" :  for  1=1  to  o;for  .7=1  tc  15;J+15I-15-K 
83:  ntf  (A5['),4'<+l,4.<+4] )  *pf.3,I,Jl  ;i"oxt  J;noxt  I 
00:  G-*X;  for  1=1  to  15 

91:  ctr  (■\S[0,3ol+4,<,3G4r4Ki)^i[I|;l+::-^;if  I>5;4to  *UoXrM 
02:  stf  (^[d,351+4<,354H;f|)*lf[0,I]  ;l+X-*X;if  I>3;3to  "uc;<i" 

03:  itl  (,VSfO,35H- W,3G?M<|) -a[2I-1| 

04 ;  ici  ( '3 [:'.,3'3++\ ,  3  jht-L I  (21) ;  1+K*'\ 

03:  "  .’eX  l":ne>.t  r 
9 G :  stx  (.'$ [0,- o7 ,4uu) )  -*-G[G  I 
07:  ret 

03:  "C;  v"  :fcr  3=1  tc  3? cot  i=l  tc  3;. •05(6(0-1.) i  3) -l..)5*X;l>V*2 

33:  for  r~l  tc  3>1 ; iX+Gfl )■••/:  GitdfLI •*■<.  nioxt  1 
.100:  tn^Z-tod-v.;  it  Z<=>;0»G  (7,  )| ;  i.cr  2 
101:  lca(:’.)-0[  3,71 
102:  next  .i;ne/.c  o 

iC3:  .vd3  .  ,13,"~itted  Jurve" ;wtc  . i,l0,D;foc  J=1  tc  6; for  1=1  tc  5 
id:  if  7(r,d!--0;./tc  ",12,”  "rj.rn  2 

105:  ore  "  !.C",efI,J] 

ICO:  next  7;. /to  i,10,13;rext  J 
107:  next  v 
103:  wtn  .’  ,13 
300:  etc  "  ' r* i •["’ 

110:  end 
*  20120 
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7.  Integration  of  dN/dr  over  r  from  0.2  to  7  microns 


This  program  integrates  dN/dr  over  r  from  0.2,  0.3,  0.4,  0.6, 
0.8  and  1  to  7  microns  for  any  specified  file  on  the  data  tape.  The 
integration  employs  the  polynomial  fit  and  the  trapezoidal  rule.  The 
matrices  serve  the  same  function  as  in  previous  programs. 

A  program  listing  follows: 

0:  "TNYEERAKS  O/VR  d’/cT  0.2  TO  7  'UtP.013  GIVI  S  TOTAL": 

1:  "A  14  VARIOUS  RAO  1 1  INTERVALS.  FROM  DA'IA  O'l  FILE" : 

2:  Jim  A$[*60| ,ES(36] ,U[6] ,3$(10)| 

3:  din  Y[6] rX[6) ,T(6,151  ,H(lol ,0(6,15] 

4:  dim  B|8| ,0(6,15] ,0,C(3] 

5:  din  R(6,16| ,Af 10) ,F[10,10] ,Mf G1 ,G[10] 

6:  715-M;dev 
7:  7-C 

3:  C+1-I:r>n  *(11,0(11  ,F[I,I] 

9:  d3D  "Insert  data  tone, continue" ;stn 

10:  ent  "Taco  *",C(3];ont  "Track  J 

11:  dsn  "fr  inter  cn,  set  form,ccntinue";sto 

12:  wtfc  M,27, 69,27,34 ,32,32,32,32,27,77,27,76,15,0,14 

L3:  v/tb  r',27,79,4  ,43 ,6,32 

14:  rnt  0,"TAPE  *",£2.0,"  INlECmTICK  OVER  RADII." 

15:  wrt  "M.0" ,C(31  ;fnt  0 

15:  wet  T:.0","F5I4  IS  FR1-’  DEDICATED  RADIUS  TO  7  MICROS  .ORTIS  .ARE  Ct-f-3" 
17:  "  Date  Tine  File  5um(.2)  Fun(.3)  3m(.4)  "-0$ 

13:  i’Vr‘i\r, :(.>'»  3 ir(.d)  Sur (l)"*i'$;wtD  .-.fL'Rf10,13 

19:  3*X  |1|  ♦-,[  21  -Xl  3]  ;1-\(51  -X  fft|  ;15-X|4| 

20:  .0345*t'fl,l] ;  .Or? 75*9(1,2] ;. 0905-9(1, 31  ;.094*P(  1,4 !;  .098-EU,5] 

21:  .102-Pf  1,5|:.  1065*Rfi,7|  :  .lll-hi|l,8] ;  .1155-Rl  1,9]  12-^(1,101 

22:  for  J=U  to  15;R(1 ,d-ll+.005-R[l,J] ;next  J 
23:  for  J=1  to  15;(.23+.025(J-l))/2-R(2,J] 

24:  .  2+.02(J-l)-R(3,J] ; .3+.00(J-l)*R(4,J) :.25t.25(J-l)-R(5,J] 

25:  if  R(5J!>l;.767*R(5,Jl+.233-F(5,Jl 

26:  1HJ-1)  -R(6,J| ;  .767*R(6,J}+. 233-R(6,J];next  J 

27:  for  1=1  to6;for  J=1  to  15;  (R[I ,J+1]+R[I  ,J] )/2-E(I ,J1 ;next  J;next  I 

28:  "start":ent  "Oeqin  with  file",E;ent  "Erd  with  file",F 

29:  for  D=£  to  F;trk  W;ldf  D,A$,E$;cll  'sort' 

30:  for  1=1  to  6 

31:  for  J=1  to  15;T(I»JI-OfI,J];next  J:rext  I 

32:  for  1=1  to  3;for  J=1  to  3;0-0(I,J| ;rcxt  T;rext  I 

33:  0*0(5, l]-3(6, 1] ;for  J=1  to  15 ; 0-0(4, J] ;next  J 

34:  ’‘nat":irc  A,[’;for  1=1  to  6;0-rQ*rl*r2-r3*rl*t5-r3-r7-r3-r9 

35:  for  1=1  to  I5;l-3;if  1=6  and  J#l;l-r6 
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36: 
37: 
33: 
39: 
40: 
41: 
42: 
43: 
44: 
45: 
46: 
47: 
43: 
49: 
50: 
51: 
52: 
53: 
54 : 
55: 
55: 
57: 
53: 
59: 
50: 
61: 
62: 
63: 
64: 
65 
56: 
57: 
63: 
69: 
70 
7L 

72 

73 

74 

75 
75 
77 
73 
70 
90 

•:1: 

>-> 

'll 

15: 

15 

05: 


if  1<15;  jira  5 

if  T>1  aid  I<6nto  "nCXt" 

if  I=l;-1.5-»X;qto  "1" 

if  1=6  and  COM  and  r7#0;Log(r8/r9)+X;log(r7/r9)*P;l»o;gt£:  "3" 
1.0X:qto  "1" 

if  OfI,Jl=0  and  r5=0;gto  "ntXt" 

if  r3t9;  ato  "cxtrac” 

if  O[I,.l]=0  and  r6=l;qto  "extrsp" 

log  (Cfl ,J!  )-F;loo(f.fT  ,J]  )*X 

if  1=1  ano  J<0;gtc  "3" 

if  I*6;qto  "2" 

gto  "3" 

"cxtroo"  :r'>l+c9;  r7+3  [1,0]  ;r  J+O  [  I ,  j  ]-*-rl;goO  "nEXt" 

"3" : r0 +>C<-rO ;  r  1+X -ri ;  r2+l-r2 ;  r 3 +PX-r3 ; r  4+Or  4 ;  if  r5=0;gto  ,,2,, 

X-9[7]  ;P*r:[ e] ;  1.5-X 

"  1" :  (r  3t2-r 4rl )  /  ( r0r2- rlrl )  +  i  [  2 1 

(rOr4-r3rl)/(rOr2-rlrl)  +  ![li;Hl]+A*:J[2]-P 

"2": for  .<=1  to  C;C+l-K-?qG+F[?,?]-F[?,n] 

W.[R]+AfP]  ;t**X-P;3X-*G 
if  K4C;d+F[R-l,R]-*Ff^-l,I3] 

3*X-*G;rext  i( 

"n£Xt":next  J;next  I 

for  1=1  to  C+l 

for  K=1  to  int((C+l)/2) 

if  N1  anc  I+XCC+2  and  I-K>0; F[ I ,1 1 -F{ I-K, I+K]  -F[ I+X  ,1  -K] 
if  I+K«+l  and  I-K+1>0;  F  ( 1 , 1+1 ] -F [ I+S ,  I-K+l  ] 
if  I+K<=C  and  I-iO0;F[I,I+l]-Ff  I-X,I+K+1) 
next  X;next  I;inv  F*F;mat  FOG 
qto  "zisioit" 

"oort":for  1=1  to  6;for  J=1  tc  15;J+15I-15-K\ 
stf  (A5 f 4X4-1, 4X+4] )  -*-r [I , J]  ;rcxt  J ;next  1 
OK;  lor  1=1  to  16 

stf  (7 5 f  261+4-< , 364+4X] ) -HI] :  1+X-K; if  ]>5;  ]to  "loXT" 

stf  (AS ( 361+ K , 354+4K] )  -n'  f  I  ] ;  l.+K-*6;  if  I>3;gto  "teXT" 

it t  (A$ [ 36]+« ,352+45] )  -a [21-1  ] 

itt  (-W f  3C3+4K  ,364+4.\J )  —X[2I] ;  1 +K—X 

"  Jc<i":roxt  I 

stf  (7$  [457 ,460| )  -r 

ret 

"zieis  it" : 

ina  J;tor  f^.2  to  6.99  bv  . 01 ; 07 -»Z ; loa  ( K ) -A  ;log  ( F+ .  Cl )  ♦  c 
for  K=l  tc  0-1 ;  Y^+5  [  K  ]  -i;  23+G  [K |  -Z ;  next  X 
t r.“Y  *'t ;  tn~<>2; .  00 5*  ( Y+Z)  4+u  [  1  ]  *U  [  1  ] 

if  v>.2C;AHJf21>ur2l;i£  •-'>.39;/'+U[31 *u[3]  ;if  F>.59;7.-HJ  [4]*U[41 
if  F>. 79; ' +G  f  51-Uf 5] ; if  1?>.99;A+U[61*i;[6] 

Text  F;fnt  G,f2.3,V",f2.9,"/",£2.0,f4.0,":",f2.0,z 
Cat  l,f5.G,fll.2,fll.2,fll. 2,£L1.2,fll. 2,  £13.2 
'Tt  ,Y [ 2| ,'/[  3 1  ,Y [1 1  ,Y [ 4 1  »x  f 5] 

vrt  -••.i«fofL|ij,.Jj2|,ui31,U|M  ,Jf5]  ,J  15] 
next  n;  jto  "start" 

'vi 


*20370 
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8.  Integration  of  dN/dr  over  r  from  1.5  to  7  microns 


This  program  is  identical  to  the  previous  program  except  that  it 
provides  integrations  of  dN/dr  over  the  intervals  of  r  from  1.5,  2,  2.5,  3, 
3.5  and  4  to  7  microns. 

A  program  listing  follows: 

0:  "l'AS-32:  imQ3,\TES  i-TO M  1.5,2, 2. 5, 3,3. 5,4  TO  7  UCP04S": 

1:  di m  «V$f 4601, C5[ 36], Uf6 1,3$ [1001 
2:  ciT  Yf6),X(6),T|S,15)  ,11(15), 0[S,15] 

3:  dim  0(M1 ,015,15]  ,LV-f31 

4:  Jim  rf5,15l,\fl0]  , F f  10, 10]  ,-’f?>l  ,G{10) 

5:  715*  !;4ev 
6:  7*C 

7:  C+l*I;rm  Mll,C[U  ,F(I,I] 

3:  dm  "Irtrorf  data  taro,ccntinus";stc 
0:  cnt  "I'arc  ,C| 31  ;cnt  "Trad;  (?",.< 

10:  Mm  "Ftinter  on,  set  forM, centime" ;rtn 

LI:  wtb  ‘,?7,5£, 27, 8 1,32, 32, 32, 32 ,2 7, 77 ,27, 76, 15,1, 14 

12:  wtb  .'",27,70,4,40,5,32 

13:  fmt  0 , \?L  r",f2.0 ,"  L'TST^ATia,  0\^r  PiiMI." 

14:  wrt  u;l , C f 31 ; trot  0 

15:  wrt  "r.3","GU:i  IS  FPJi  LTDICArED  FADICG  1C  7  ;  iICEX  15  .UNITS  APE  C^-3” 
16:  "  Date  Time  File  Sirod.5)  Sur.(2)  3uti(2.5) 

17:  3$.V’Sur(3)  Sum(3.5)  Sun(4) ”-3$;wtb  10,13 

18:  3*X[11*:<(2]*X[31  ;1*X(51*X[61  ;15-»X[4] 

10 :  . C‘345+2 fl,ll ;  .0875*4(1,21;  .000 5*  l [  1,3];. 0S4*R (1,41;. 013*R (1,3] 

20:  .102*r;fl,G|;.lC55*'7fi,71;.lll*P(l,3|;.1153*ra,l/|  ;.12+*11,1 0) 

21:  for  3=11  tc  i6;Ofl,J-i)+.003*:;[l,J]  ;ncxt  J 

22:  ror  3=1  to  15;<.23+.')23(J-1))/2*7[2,J1 

23:  ,2+.u2(J-l) *^( 3,Jl ;  . J+.06(J-i)*P(4,JJ ;  .25+.25(J-l)*F(5,J] 

21:  if  3(5,3 1 >1; . 7t>7*P(  5,3] +•  233*7(5,3] 

25:  l+(  I-1)*3(6,J1 ;. 767*0(6^1]+. 233-*7(5,Jl  ;next  J 

26:  for  1=1  to  6;  for  J=1  tc  15;  (R(I,.T+lj+f:(I,Jl  )/2*fJ(I,  J];rext  J;r.ext  I 
27:  "rtort"  :ent  "fxKjin  witn  £ile",Z;enc  "Tru  witn  f  ile" ,  F 
28:  for  n=F.  to  r  :trk  .v; lof  D,A$,C$;cll  'sort' 

23:  for  1=1  tc  5 

30:  tor  1=1  tc  15;T[I,Jl*0(I,J];rext  J;next  I 

31:  for  1=1  to  3;  for  .7=1  to  3;OOfl  ,J]  ;next  J;next  I 

32:  00(5,11-0(5,1]  ;for  J=1  to  15;0*O]4,J]  ;next  J 

33:  "Mat "ana  A,F;for  1=1  to  5:0*r0*rl*r2*r3*r4*r.5*r6*r7*r3*r9 

34:  ter  J=1  to  15; 1*3; if  1^5  ani  Jrl;l*rb 

35:  if  J<15; ior  5 

35:  if  [>1  me.  I<6;qto  "nrXt" 

37:  if  I=l;-1.5*X;ato  "1" 
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] 
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38:  if  1=5  md  rO^O  and  r7  *  0 ;  lcc  ( nJ/cJ )  +X ;lcxj  <r  7/r9 )  -F;  1- r5 ; -3 to  "3" 
39:  1.5*X;gto  "1" 

10:  if  O[I,7)=0  crvi  r6=0;gto  "nh'Xt:" 

41:  if  r9-iO;  Jto  "ex tree" 

42:  if  OfI,Jl=0  ana  r6=l;gto  "cxtrao" 

43:  l»(OfI,Jl  )-F;log(r[I,J)  )+X 
44:  if  L=1  ana  a<9; gto  “3*' 

45:  if  liojcto  "2" 

46*  c'to  11 

47:  "extr so" : c 9+  l*tcj ;  r7  +3  f  I ,  J  7 ;  r  6+3  [  I ,  J  l-*-r3 ;  alo  ,,nrXt“ 

43:  "3"  :r0f.s;<+n');rl+»rl;r2+l*c2;r3FFX*r3;r4+l-*r4;if  r5=0;cto  "2" 

49:  »Lf7]  ;P*fi[(jj;1.5*X 

50:  "1" :  (r3r2-r4rl)/( r0c2-rlcl) ■*  1[2) 

51 :  ( r0r4-r3 rl )  /( rf)r2-rl rl)  -rf-tf  1] ,  m  [1]  +X*Mf 2 1 -F 
52:  "2" :for  K=0  to  C;C+l-K-*«;G+F{P,Rj*F[R,Rl 
53:  P+AfF|*AfP];P*X+P;GX-*G 
54:  if  Kf£;J+FlR-l,R|+F(R-l,R] 

55:  G*»G ;  next  K 

56:  '‘nCXt":next  J;next  I 

57:  for  1=1  to  C+l 

58:  for  X=1  to  int{(C+l)/2) 

59:  if  1^1  and  l+KCC+2  arc!  I- 4>0;f[T,I j ♦rfl—X, I+4]-»i* [ I+ICr I — .-») 
aO:  if  I+K<=C*-1  anJ  I-i;+l>0;FfI,IH]*FlH-K,I-.-;+l) 

61:  if  I+K«  anJ  I-K>0;F(I,I+lJ*tfI-:\,I+K+ll 
62:  next  K;ne:<t  I;inv  F+Fjirnt  FA*G 
83:  ate  "eisisit" 

54:  •’eort,‘:£er  1=1  to  5;  for  J=i  to  15;  J+15I-i:>*.< 

65:  r.tf  (•’'514*’+1, 4K+4))-*fTlrJ)  :r,ext  J;next  1 
56:  0 ->■■■;  for  1=1  to  15 

67:  stf  P\?r.r>lt«,354f4,v| ) -Hfl) ;1+K+K;if  l>5;atc  "ic.'CJ" 

53:  ett  (4 $ f  351+ 4X , 364 +K1  W  ( 1 1 ;  l+:>5; i f  I>3;ito  •‘teXT” 

69;  itf  f  3«1  +4<  ,352+4 <1  )-<f2f-H 
79 :  itf  ('0 f  35 3+4X  ,364+4X1 )  ♦«<( 211  ;l+\’-»i\ 

71:  ”'le'(T":no.xt  I 
72:  atl(AS(1'7,;40|)-f 
73:  ret 
74:  "zisieit": 

75:  iaa  d;for  P=1.5  to  5.9J  hy  .02;0-»Y*2;loa (?)*A;l0'!(F+. 02)*2 
75:  tor  4=1  to  C+l;Y/'+5[<(lvi';7,B+j(ts|*z;rext  K 
77 :  tn/‘f-*Y ;ta~L+Zj  .01*  (Y +4 )-»A;A+J  fl|  V)  [1] 

76:  if  r>l. A+Uf 2] •*u[21 ; if  R>2.48;A+U[3]-U[31  ;if  R>2. Oci; A-H j f 4 ] -MJ f 4 1 
70:  it  P>3 . 4(3 ; A+u [  5]  -HJ  f  5 1 ;  i  f  F>3.?3;A+U[6]-o[C| 

80:  next  F;f.ot  0,f2.0,"/%f2.0/VV2.9,f4.0,  V',£2.0,z 
81:  Lot  l,to.c,112.4,fll.'4,fU.4,ill.4,fll.4,tll.4 
»2:  wrt  ,Y( J|  rx  [1]  ,f(4|  [ 3 1 

33:  wrt  1. 1'  ,.,ufl] , J 12) ,6[3| [4|  ,L’[5|  ,J  [5| 

84:  next  I;;  j  to  "rtart" 

35:  eru 
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9.  Plot  Fractional  Change  in  dN/dr 


This  program  plots  (l/dN^HdNg  -  dN^)  as  a  function  of  log(r) 
from  the  data  in  any  two  specified  files.  It  does  so  by  utilizing  poly¬ 
nomial  fits,  in  log  space,  to  the  two  data  sets.  The  matrices  employed 
have  the  same  correspondence  as  in  previous  programs. 

A  program  listing  follows: 


0:  ’VA2-32:  ADVZ  Fn/CilOJAL  CH/YKTS  IN  dl/dr. ..  (14.*)  (>3r.:LlA fili/dr) )": 

1:  Jim  '.$[2, -’.301 ,0$f4G01  ,FS(36]  ,  J  (21  ,ofo,15(  ,Cv(25|  ,Z[2] 

2:  dir  Y[2,31  ,Xf61,rf2,.3f15),;)|I3|,0[2,o,151 
3:  Jin  .:(31,:-:(6,151,J,CI3],X[10| 

4:  din  OfS,15|  ,9(3,161  ,4(101  ,Ffl0,101  £[101 ,1  f91  ,M[6I 
5:  35*3(1| ;11I*3[21 ;43*d(3)  ;Q*6f4|  ;  42*8(51  ; 64*3(6) 

6:  715 *4;dcv  "T,M 
7 :  fnt  0,10x,z 
B:  fmt  1,0.0 
9:  int  2,cl,f5.1,z 

10:  fmt  3,"Cate",f3.0, ,f2.0,"/",f2.0, "  Time",f3.0,":",f2.C,z 
11:  fmt  4 /'Averaging  tine  *  ",f2.0,"  and  ",£2.0,"  minutes" 

12:  fmt  6,el0.2,z 

13:  fmt  7,"Tsne  #",f3.0,"  Piles”, £3.0 r n  and",f3.0, "  Fvcnt  ",c25 
14:  fmt  8, "Polynomial  of  order  ",f2.0,z 
15:  fmt  9,el5.7,z 
16:  7*C 

17:  C+1-I;rdm  Afl] ,3(1) ,K(I]  ,P[I,I] 

18:  dsn  "Insert  data  taoe,cantinue";sto 
19:  ent  "Taos  #",C(31 ;ent  "Track",,# 

20:  dsp  "It inter  cn,  set  form,ccntinue'';stP 

21:  wtb  M,27, 69, 27,34,32, 32, 32, 32, 27, 77 ,27, 76,15,0 ,14 

22:  wtt>  J-',27,79,4,43,6,32 

23:  3*X[11*X(21*X(31;1*X(51*X(61;15*X[41 

24:  .0845+R[  1,11 ;  ,0375*Rfl,21 ;  .0905*R(1,3) ;  .094*R(1,4) ;. 098*R[1,5) 

25:  .102*8(1,61 ;. 1065*R(1, 7]  ;.m*P(l,8);.n55*R(l,9J;.12*Rll,10) 

26:  for  J=ll  to  15;R(1,J-1)+.005*R[1,J] ;noxt  J 

27:  for  J=1  to  16; (.23t.025(J-l))/2*R(2r3l 

28:  .2+.02(J-l)*R(3,Jl ;  .3+.03(J-l)*R[4,J| ;  .25+.25(J-l)*Rf  5,J] 

29:  if  R(5,J|>1;  .767*R(5,Jlt.233*R(5,J) 

30:  l+(J-l)*Rf5,Jl  ;.767*P(6^1+. 233*1^(6, J1  jnext  0 

31:  for  1-1  to  5: for  J=1  to  15; (B(I,J+li+8[I,Jl )/2*E(l,JJ;mxt  Jjnext  I 

32:  ''rtart":ont  "f Vent", CS; ent  "'"ackqround  file",Jfll  ;cnt  "LVcnt  file",J[2) 

33:  uck  7;for  7=1  to  2;llf  J(i7l , 3S,CB;2.$*A3(L'I ;cll  'scrt';rcxt  C 

31:  for  to  2; Lor  J=L  tc  6 

}5:  l or  J=1  Lo  15; f(L,I ,J|*C(n,I,Jl ;rext  J;noxt  I 
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35:  for  1=1  to  3; tor  J=1  to  3;C-*0[C,t,J];next  J;next  I 
37:  0-*O[0, 5, 11-0(0,5,11  ;fcr  J=1  to  l3;0-*O[O,4,Jj;next  J;;cxt  C 
35:  for  1=1  to  5: for  J=l  to  15;.}(2,I,JK)(1,I, JJ 
39:  if  :'(I,J  ]<U;OS(I,J| 

-10:  next  J;next  I;for  i>l  to  2 

41:  "mat":ina  '\,F;fcr  1=1  to  5 ;  O-OJ ♦rl-*-r2*r3  +r4-*r5-»r5 -»r7*r 3*rQ 

42:  for  J=1  to  16:1-0;  if  1=6  ani  J#l;l*r6 

43:  if  J<15;yr.o  5 

44:  if  I>1  anc  I<5;qto  "nEXt" 

45:  if  I=l;-1.5*<;gto  "l" 

46:  if  1=6  and  r9i(0  ard  r71G;  log  (r6/r9)-*X;loa  (r7/r9)-*-P;l-*-r5;gto  ”2' 
47:  1. 5*X;  gto  "1" 

18:  if  0(D,I  ,J]=0  art!  r6=0;gto  "nEXt" 

19:  if  r9#0;gto  "extrao" 

50:  if  0[D,I,J)=0  and  r6=l;gto  "extrao" 

51:  log  (Of D ,1 , J|  )-*-P;log  (F. [  I  ,J] )  ♦X 
52:  if  1=1  ana  JO; gto  "3" 

53:  if  I  iff;  gto  "2" 

■54:  gto  "3" 

55:  "ex tr oc " :r 5H- l-*r 3; r7+-0[J,I,JJ-»r7;rt?+E(I,J]-»-r3;gtc  "nEXt" 

55 :  "2' :r0Ev<-*rJ;rl+X-»rl;r2+l>r2;r3+PX*r3;r4+P*r4;if  r 5=0; gto  "2" 

57 :  X-O  [71  ;f-*  '  f 6| ;  1. 5*X 

56 :  " .1 " :  ( r 3 rl- c4rl )  /  ( rO r 2- rlrl )  -*■  i[  2 ) 

53:  ( rOr4-r3 rl )  /( rl)r2-rl rl) -«' '( 1]  ;"(ll  fX*i:(2|-0 

60;  "2"; for  1=0  to  C;  2-il-h-*-*;  GO 

61:  P+a  ft  1  r rl ;  P*X»P;  2<-*C 

52:  if  K=c:;3-H-(R-l,3K-fP-l,Pl 

63:  2* no xc  K 

64:  "nEXt" :next  J;next  I 

55:  for  1=1  to  Ol 

65:  for  '=1  to  ir.t(  (C-rl)/2) 

57:  if  I,fl  and  XO<3+2  and  l-0>0;r'f  I,I]*fT(I-i<,I-H\)-*r  (l+K,i-l] 

66:  u  I+t«C+l  ard  I-i\<-l>0;P(I,I+l] -*■!*[ I+K,I-K+11 
69:  if  I+K«  and  I-iOC;F(I,I+ll*r’[I-K,I-t-K+l] 

70:  next  F>’:next  I:inv  F*F;  if  D=l;rat  FA*G;iaia  2 
71:  net  FA*K 
72:  next  D 

73:  "out":wrt  "m.7",C(3]  ,J[2]  ,Jfl]  ,C$ 

74:  Y[2,4]-Y[l,4]*Z;Y[2,5j-Y(l,5]+60*Z*Z 

75:  wrt  "M.3",Y [2,2]  ,Y(2,31  ,Y[2,1J  ,Y[2,4)  ,Y(2,5] 

76:  fmt  5,"  u'i-",f3.G,"  minutes"  ;wrt  "M.5",Z 

77:  wrt  'H4",Z(2»,ZflJ ?wtb  M,10,13 

78:  wtb  iv,27, 55,-4 ,0,7,32,"  (1/1)  (!3ELTA(di}/ciR) ) " 

79:  wto  '1,27,65,-4,-210,-1,-16, "log(radius) " 

30:  "plt":if  flg7;gto  "skin  orint" 
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31:  -1*X; -4*Y 

82:  wtb  19,27,65, irt(15X/4) ,int(24QX) ,0,0 

83:  if  X*0  and  XnxxH=0;wtb  M,"|  ",10,3 ,3,8 ;wrt  "M.l",X;qto  +2 

84:  wtb  M,"-" 

85:  if  (X+.05*X)  <2.0$;4to  -3 

36:  wtb  m,27, 65,0,0, int(1.5Y) ,int(9SY) 

87:  if  WO  and  i’n.ccl=0;wrt  "M.2","-",2*Y;-jto  ¥1 
83:  wtb  r,"|" 

SO:  if  (Y+.1»Y)  <5.1;'itc  -3 
30:  Jto  "skin  nr  ir  e" 

91:  "3c.rt‘':for  1=1  tc  Gjfnr  J=1  tc  15;J+15I-1vv*-a 
92:  otf  (A5fC,4'C+l,1Kt-41)-*rf3,I  ,J1  ;rcxt  Jjnexc  I 
33:  0*.\;  for  1=1  to  1G 

94:  .~tf  (ASfO,3ClwK,354t-4:<l)^ifI|;lK>Kjif  r>G:?to  " feXf" 

95:  stf  (r'.$[d,3Sl+-4\,3Gi+4K])*Yl  ),Il  :l-i\>u;if  I>3ntc  ".leXf" 

%:  it!  (1$ (0 ,3C1+K' , 352+4X1 ) *< {21—1 1 

97 :  itf  (Av  !  0 , 35 3  *-4X , 354+4X 1 )  ♦Xf21 1  :1+K*K 

09:  "  5e<r’:n3Xt  I 

>9:  ntf  0810,437,4601  )-2[C] 

100:  rat 

LO 1 :  "skin  nr  int" : -1-*X 
102:  “crr:l+Y;>Z 

LOO:  ear  1=1  to  C+-1:  YX+Ctll  -*i;3XT.s(I|-»^;ncxc  I 
104:  trfd/trTx-l-*-!: 

105:  if  Y<~3  or  Y>10;  jr.o  4 
106:  •5*(i’t’Vfi6)  *i 

107:  vtn  '',27 ,63, int(15X/4)  ,ir.t(240X)  ,int(3Y/2)  ,int(9oY) 

103:  wto 

103:  "nclp":ir  (X+l/120+X)<.33;jto  "crv" 

110:  v/tn  3,12,13 
111:  ?eo  "start" 

112:  enO 
*17154 
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10.  Plot  Normalized  Fractional  Change  in  dN/dr 


This  program  is  identical  to  that  just  preceding  except  that  the 
plot  is  normalized  by  the  peak  value  of  the  fractional  change  in  dN/dr. 

This  is  accomplished  by  determining  the  maximum  value  of  the  fractional 
change  in  dN/dr  in  the  radius  interval  of  interest  and  plotting  the 
fractional  change  function  divided  by  this  maximum  value.  All  such  curves, 
then,  have  a  maximum  of  one.  The  matrices  are  the  same  as  in  the  previous 
program. 


A  program  listing  follows: 


0:  "|>\>32:  PTOiS  (1/4)  FXTTMd.ydr))  i-XK  AL12T.0  P*  ITS  PLv't.'  VAUJi?": 
1:  dir  ASf2,46Gl  ,”.•?( 466 1  3c  1  ,J  121 ,3(0,151  ,C5(2 5]  ,1121 

2:  Jiir  Y(2/.l  r(2, 5,15), :M15],;(2,-5,15] 

3:  lin  i(o), 2(6,151 ,n,2(3l, 4(10] 

4:  diu  1 6 , 1G )  ,!>[o,1j|  ,\(10|  ,r  (10,101  ,JflO]  ,1 (9]  ,  ,(61 
5:  35*7(11  ;lll-u-(2l  ;43-il31  ;^5(4]  ;42-:;(5!  ;64*u(6) 

5:  715-  t;dcv  “  * 

7:  Llit  0,lCx,z 
5:  uit  l,L3.d 


3:  c . t  t~ ,cl ,  - j •  1 , z 

10:  nrt  3,"5.ite",f3.j,"/",f2.0,"/",i:2.0,"  Tire" 
11:  fnt  4 /'"vcrcniPQ  tire  =  ",12.0,"  and  ",f2.0 


,f3.0 ,":",£2.0  ,z 
,"  minutes" 


12:  tmt  5,elQ. 2,z 

13:  fmt  7,"Tace  S“,t3.0,n  Files ",f 3.0,"  arc", f 3.0,"  tvent  ",c25 
14:  ent  3,"lolvnorrial  of  order  ",f2.J,z 
15:  fait  9,el5.7,z 
16:  7-C 

17:  C+1-I;rdm  AfI),G(I] ,K(I] ,F[ I ,1] 

18:  dco  "Insert  data  tare, continue" ;sto 

19:  ent  "Tane  #",C(3];ent  "Track", W 

20:  ;1sd  "Printer  on,  set  form, centime"  ;stn 

21:  wtb  f*,27, 69 ,27, 84, 32, 32, 32, 32 ,27, 77 ,27 ,75 ,15,0 ,14 

22:  wtb  f',27,79,4,43,6,32 

23:  3*X(ll*Xf2I-Xf31;l*X[5]-X(5];15-tf[41 

24 :  .  0345*R  f  1 , 1 J ; .  087  5*R  [  1 ,2] ; .  0905*R  [  1,3 ) ; .  094-R  [1,41  ; .  098-R  [  1 ,51 
25:  .102*r(l,6J;.1055*R(l,7);.m*P(l,3(;.1155*R{l,91;.12*R[l,101 
26:  for  J=ll  to  15;P[l,J-l]+.005*R(l,Jl  ;next  J 
27:  Cor  J=1  to  15;(.2>.025(J-1))/2*«(2,J1 
28:  .2+.02(J-l)-R[3,Jl;.3+.08(J-l)-R[4„J];.25+.25(J'l)-n(5,J] 

29:  if  P(5,J!>1;  .767*R(5,Jl+.233*Rf5,J) 

30:  it(J-l)*f’f6,J)  ;.757*R(S,J]+.  233-*R(6,  J)  ;next  J 

31:  tor  1=1  to  5;for  J=1  to  15;  (R[I,j+l)+R[IrJl  )/2*2(I,J(;next  J;naxt  I 

32:  "start" :ent  "f vent" ,23; ent  ".'ock.jrcuru  file"  ,J[11  ;ont  "tvent  rile",j[2] 
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33:  trk  '.J;fcr  !>1  to  2;ldf  J[D)  ,3$,i3$;'33^(Dl  ;cll  'sort'jrcxt  0 

34:  for  D=1  to  2; tor  1=1  to  6 

35:  Lor  >1  to  15;T[D,I  ,Jl-*C[D,I,J|  ;next  J ;noxt  I 

36:  tor  1=1  tc  3;  tor  J=L  to  3;0-*Ofr,I,J]  ;noxt  Jjnext  I 

37:  00 f  0 ,3 ,  U  -»C  p  ,5 , 1 1 ;  f  or  J=1  to  15;Q*0|i:,4,)];next  J;rmxt  0 

33:  tor  1=1  to  6; for  J=1  to  15:Jf2,I,J]-0[l,I,J |-3 [£,J] 

33:  if  SlI,Jl<d;i**S[I,J| 

40:  next  J;ncxt  I; for  ;>1  to  2 

41:  "rat"  :ina  A,F; for  1=1  tc  5;G-*-rJ-*rl-»r2-*t3 *r4-*’r5-»-r6-*-r7*rG+r3 

42:  for  J=1  to  16;  1*44; if  1=6  orvi  J*l;l-*rS 

43:  if  J<15;ioc  5 

44:  if  I>1  and  I <5;: to  “next" 

45:  if  I=l;-1.5*:intc  "l" 

•"•6:  if  1=6  ana  r9r0  arri  r7 40;loa (r3/r9) ♦X;lon(r7/t9)  +P; l-*-r5;tto  "3" 
47:  1.5*X;gto  "1" 

43:  if  Op,I,J)=J  art  rbO;qto  "r£Xt" 

49:  if  r9*0;qto  “extraD" 

59:  if  0(D,I  ,JJ=Q  art]  rb=l;gto  "extrao" 

51:  loq(0[D,I,J)  )-»P;loa(E[I  ,J]  )-*A 
52:  if  1=1  and  J<9;uto  "3" 

53:  if  I46;crto  "2" 

54:  qto  ” 3" 

55:  "oxtrao":r9tl-*r9;r7+OfO,I,J|-»r7rr3tF[I, J]*r3;gto  "next" 

56:  "3":rO+XX*rO;rl+X*rl;r2+l*r2;r3+PX*r3;r4+P*r4;if  r5=0;qto  "2" 

57:  X*3  f  7] ;  P*r  f  P] ;  1 .  5-*X 

52:  "1":  <r3r2-r4rl)/(rr.t2-rlrl)v .[21 

39 :  (rOr4-r.3 rl)/( n»r2-ri rJ )  v  flj  [ l)  ix*: >12) *;■ 

60:  "2": tor  ’C=)  to  C;C+l-5*0;C^f |F,Pl*?fn,P] 

61:  o+a 

62:  if  aFC; JtFfp-l  ,i'.|  ♦fil’-l ,3 1 
5.3:  6*  X-*3  •  next  \ 

64:  "nLlLt"  :ncxt  4;. text.  [ 

55:  for  1=1  to  C+l 

65:  for  \=L  to  mt(  (C-hU/C) 

67:  if  I  pi  and  i+K<Ct2  ana  I-0>0;L'[I,Il-*f[I-r\,Itk]-*T[l+i',I-»<] 

63:  if  H-K<=d+1  and  1— \KlX,;F[I,T+].]-*v[I+-X,I-K+il 
55:  if  I+K<=d  an.]  I-0>  0;  f  { I ,  I  t-l|  [  I-u ,  I+K+J  ] 

70:  next  A;noxt  I;inv  F+F;  it  l>l;.nat  ?A-*d ; y,n  2 
71:  ~et  i7*<\ 


72:  next  .3 


73:  "cut":wrt  ,M.7",C[3|  ,0)21  ,J|  11,35 

74:  Y[2,4)-Y  11,41  *2;Y  (2,5]-Y[l,3]+6G*>Z 

75:  wrt  i.3”,Y 1 2 ,2)  ,Y ( 2 ,3]  ,Y|2,1J  ,Y[2,H  ,Y [2,51 


76:  fart  5f"  ll=",f3.9,"  ninuter" ;v r t  U'.5",Z 
77:  art  f21 ,3(11  ;./tb  ’,10,13 

T-'i  'Jtb  '  ,2  7 ,55  ,-4 ,0 ,7, 32 , "<  (1/  In)  (l/i)  p;  tOMd/oP))  >” 
70;  atr>  27, 65,-4 ,-210,-1,-15, "lcxi(ro.lit*s)  " 
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30: 

"rlt":i 

t  f 3  :7;  }to  "ok i~>  or ir.t" 

31: 

-»v 

82: 

7,09 

,i;it  (]  3.v/‘!)  ,  int(24CX)  ,0,3 

S3: 

ic 

OHi; 

Krail=0;wf>  V | ",10,8,3,3 

;-.-Tt  "'*.1",K; qta  +2 

3-1: 

Vtt- 

^•1 

S3: 

it  (X+. 

OW 

)  <2.05;qto  -3 

63: 

wtb  :  >2 

7 ,5  j 

,0,0,int( 1.5Y) ,int(96Y) 

97: 

it  YrO 

nnw 

i::oil=0;  ;rt  "  1.2", 

Mto  +2 

63: 

.•/tb 

i" 

£9: 

it  (Y+. 

NY)  0.1;  7 to  -3 

90:  itc  "nkio  nrint" 

SI:  "3crt":for  1=1  to  3;£or  J=1  to  13;0ti5I-16*X 
92:  st£^S[0,4;<+l,4i<+41)+p[0,I,JJ  ;ncxt  J;ne>±  I 
93:  OK;  for  1=1  to  16 

94:  stf  (AS [") , 361+4K , 3j'4+4,<] )  -*0[  I ) ; ItOK; if  I>5;nto  ’lieitf" 
95:  r.tf  (AS|0,361+4X,364+4K])*Y[0,I]  ;l+K*K;if  l>3;qto  "^eXT" 
95:  itf (AS [ 0 , 361+4K , 3o2+4K ] ) -X [ 2 r-1 ] 

97 :  i tt  (AS fD,353+4K , 364+4KJ ) -X (21 J ;  1-H<*X 

98:  "A'cXT'^next  I 

99:  3tr (AS[0,457,460))*Z(01 

100:  rot 

101:  "skin  ori:it":9*A;ror  ^=1  to  2;-l*X 
102:  "crv" 

103:  £cr  1=1  to  C+l;  yX+J(I|-Y;3K+;<H]*Z;next  I 
104:  trfZ/tn~Y-l*i ;ir  9=2;~nm  3 
105:  it  Y>A;x’-*7' 

106:  3  to  "helr'" 

107:  Y/.VY ;  it  Y<-1  or  Y>l.i;ir..n  4 
103:  5*Y-Y 

109:  wtb  M,27,b5,int.(15X/4)  ,int(240X)  ,int(3Y/2)  ,ir.t(96Y) 
110:  wtb  -VV* 

111:  "bole" :i£  ( <+l/i20*;\)<..33;tto  "crv" 

112:  next  0 
113:  nth  ‘!,12,13 
114:  qtc  "stort" 

115: 

*37'4 


1 
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11 .  Tabulation  of  Aerosol  Decay  Constant 

This  program  calculates  the  decay  constant  for  aerosols  as  a 
function  of  radius.  It  assumes  an  exponential  decay  and  requires  the  input 
from  three  data  files,  the  third  being  the  background.  The  time  elapsed 
between  the  two  post-event  data  sets  must  be  specified.  The  calculations 
proceeds  from  polynomial  fits  in  log  space  and  is  given  by: 


'  i 

I 

J 

't.i 
•  ♦ 


f 

4 

f 


t  =  (t^  -  t^ )  (InCd^/dr  -  dN^/dr)  -  In  (dN^/dr  -  dNb/dr)) 

i.  L 

Where  t^  is  the  time  of  the  i—  file  and  the  subscript  b  denotes 
background. 

The  matrices  employed  are  the  same  as  in  foregoing  programs 
except  the  dimension  of  some  arrays  has  been  increased  by  one  to  accommodate 
inclusion  of  the  background  data  set. 

A  program  listing  follows: 


G:  ";PP 


iWiHATES  •i’i’K  AFJXXCL  t-CCAY  OUSCNTl',  TAU,  SOP  VARIOUS  RATII": 


1:  Jin  ASf3,.16Cl  , ’>Sf 460]  ,P$[361  ,J[41  ,C$[251  ,Z[3];sfg  14 
2:  dim  Y [3,51  ,XfG|  ,Hl  161  ,0[3,6,L51 
3:  dim  E[6, 151  ,3,C(3I  ,X[1C1  ,316,151 
4:  dim  R[6,16|  ,M10|  ,Ff  10,101  ,Cfl01  ,IflOI  ,Mf6|  ,D$f25] 


r;  • 

M 

"*D$ 

6: 

701*M;dev 

7: 

fret.  0,10x,z 

8: 

fmt  1,0.0 

9: 

fnt  2, cl ,£4.0,2 

10 

fmt  3, "Dote", 0.0 

,  £2.0,"/", 0.0,; 

11 

fmt  4  ,‘"\veragina 

time 

=  ",t2.0,”  and 

12 

fmt  6,clO. 2,z 

13 

fmt  7,  "I one  ",£2. 

0,"  . 

ikqnd  ",f3.0, "  ! 

14 

fmt  9, old. 7, z 

15 

7-*C 

16 

:  Ol-Ijrac  MU  ,3  [11 

[11  ,i[il,i(i,H 

17 

:  '‘ten  " Inner t  data 

tare 

.continue"; cto 

13 

:  ent  "Inre  # " ,C { 3  J 

;ont 

"Track",, v 

id 

v  i3D  ".-r  inter  on. 

ret 

form.caitinuc"; 

2J 

vtr>  ‘ 27, 59,27, 01 

,32, 

32,32,32,27,77, 

21 

wtb  ',27,70,1,44, 

6,32 

22 

2*'i\  1 1 1  1 21  »X[  31  ; 

1-K[ 

51  *X[6| ;15*X[4] 
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23:  .0S45*Rfl,l] ; .0875-3(1,2) ;.0905*Rf 1, 3] ;.094-K(l,4] ; .C9d-R(l,5) 

24:  .102-P(1,6|;.1065-S[1,71  ; .111-Pf 1,8] ; .1155-P( 1,91  ;.12-R(1,10] 

25:  for  J=ll  to  16;R(1 ,J-11+.005-R(1,J| ;nsxt  J 
26:  for  J=1  to  16; (.23+.025(J-l)) /2-R(2,J) 

27:  .2+.02(J-])-Rf3,J];.3+.08(J-l)-R(-l,Jl  ;.25f .25(J-1)-R[5,JJ 

28:  if  P(5Jl>l;.767*P(5,J]+.233-R(5,Jl 

29:  l+(J-l)+Pf6,J] : .767*Rf 6,J] +.233-P(5,Jl ;next  J 

30:  for  1=1  to  5; for  J=1  to  15; (^fI,J+ll+PfI,J])/2-E(I,JJ:next  J;ncxt  I 

31:  "START" : for  V=1  to  2:wtb  '-<,10,10 

32:  "rtort":ent  "Kv&->t",::$;C5SOS (l,24-3on(CS)  1  <$ 

33:  cr.t  "  TACTS:-:  (HKD  FU,3?",J(1]  :if  3(l]=J[4l;5f<4  1 
34:  J  (1) -J(4) 

35:  cnt  "2.-Cn  Fir.t,,,Jf2J:Snt  "StfU’G  FILF", J (3) ;cr»t  "Till  DIF?",r 
35:  crk  ,v;for  i>L  tc  3; if  flil;cfa  1; tto  2 
37:  iof  Jfrl ,.  J,r;$:.3S*/'Sr:->|;cll  'fCCt' 

33:  no:<t  R;fcr  [=1  tc  5:  for  J=1  to  15 

39:  2(2,1  ,Jj  -Ofl,l  ^1 1-0  (2, 1  ,J|  ;0(3,I,Jl-Ofl,I  ,J]-0[3,I,J  J 

40:  if  0(3,1, Jl=i);-lC000-3(I,Jl;qco  "BLTFAST" 

41:  0(2,I^J)/0(3,I,J)-Q[2,I,J)-3[I,JS 
42:  if  Sf I ,J  ]<=0:-]0000-S (I,J); iro  2 
43:  T/ln(P(I rJ] )*3( I,JI 

■U:  " BELFAST" :ncxt  J;ncxt  I; for  1=2  to  3;0-C(3,l,I]-O(3,2 ,1)  jnext.  I 
45:  for  1=2  to  3;-100P0-S[l, 11-6(2, I] ;next  I 
46:  for  02  to  2 

47:  "oiat":ira  4,F;for  1=1  to6;if  I=4;qto  "rFXtl" 

43:  for  J=2  to  L5;1-*G 

49:  if  2[3,I,J|=0;4to  "n(7<tJ" 

50:  0[D,IrJl-P;loa(E(T,J|)-< 

51:  "2"  :for  K=0  to  C;C+l-K*P;C-tP[P,Rl-r  (H,R| 

52:  u+A[Pl-A(Kl  J 

53:  if  2  =C;3K'(P-l,  '!-f  t2-l,M-l 

54:  C^X— 3:r.oxt  4 

55:  "nC4tJ":roxt  J 

56:  "nlXtl" :roxt  I 

57:  tor  1=1  to  C+l 

58:  for  ,\=1  tc  int((01)/2) 

5Q;  if  i si  on-i  ( fK<C+2  one  £-5>0;F  (I,  ()-(•  fI-X,I+K|-f(  I+i’,i-i\| 

Of. :  it  I+XCC+l  ?.n3  i-;oi>0;F(i,i+li-P!i+A,i-;;t-ll 
61:  if  I+.K<=C  arti  I-.<>(!;F[I.I+]|-F(I-*,I+*+ll 
62:  next  K;r.oxt  I;inv  in  t=l;rrot  FVG;>to  3 
53:  if  C=3;:;. 3 1  (VK’jr-n  2 
54:  T»t  FA- 1 
63:  next 

56:  "nuf":v,rt  ■••.7" ,2(31  ,1111 ,7(21  ,7  PI  ,f 

57:  -..rt  -  .'>",yfl,21  ,Y[1,31  ,'/(!, II 

65:  t>t  ”  ,•  v?nt",c25;wct  ”A",0?;.;to  (1,10,13 

7):  wth  '  ,"rci!  (rioutoc)  ";v/th  -’,in,l3;for  J=2  tc  15;for  1=1  to  5 


§ 
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70:  if  1=4, -ate  "CAUAY" 

71:  if  P[I, j)<-9999;wtb  V  ",42,"  ";irro  3 
72:  if  SfI,J]>9c99;wtb  M,"  ",42,"  ";jntp  2 

73:  wrt  ’M.6",S[I ,J) 

74:  "GALYA  Y" :  next  I;wtb  >1,10, 13;  next  J:wt±>  M,10,D 
75:  qto  "erv" 

76:  v;tb  M,27,55,0,0,int(1.5Y)  ,int(96Y) 

77:  if  Y#0  and  Yrcoi=0 ;wr t  ">i. 2",  "-",2Y;qto  +2 
73:  wtb  M, " I " 

79:  if  (Y+.1-»Y)  <5.1:  qto  -3 
80:  qto  "skin  crir.t" 

61:  "scrt":for  1=1  tc6;for  J=1  to  15;J+15I-10K 
32:  ntf  (‘’'$[D,4(\tl,4K+41  )-*OfC,I,J)  ;next  J;next  I 
03:  OK: for  1=1  to  15 

54:  stf  (A$ |D,351+-4K ,3f>4+4K] )-*■'!{  I ]  ;l+OK;if.  I>5;  ito 

35:  stf  (AG[:5,3r>l+-iX,354Hi\|  )+Y[0,T]  ;]fOK;if  I>3;atC  "  loXI'" 

36:  iti  (AS  [9,361-1-  ’-K  ,362+4X1 )  ♦.<[ 21-11 

37:  iff  (\$ 0,36 3+45, 36 4+4T<] )  -X[2l[  ;!+'<*< 

sO:  "i !oXi'":oc>'t  T 

W>:  stf  <  AS fO, 467, 455| )  -ZJ.Ol 

90:  ret 

91:  "skin  crir.t" : -1-*X 

02:  "erv" : t  or  ,'=l  to  3:  Cor  .7*1  to  5 ; . 05 ( 3 U-l )  +0 ) -1 . 05*vv ; O Y 

93 :  for  1=1  to  C  H :  YY-X  ft]  next  I 

04:  Y*6f.],r,i*;-.ii  /:<=,;-ir000-»-f  j  ? 

55:  Jj  j  , 4 1 

y> :  next  5 

57:  wti:  t',13, "Kitted  Curve" ;wtb  •  •,]'.), 13:  fee  J-l  co  3; for  T=1  tc 

53:  ir  5|I,.:|<-,eoq;we)  9,"  ",42,"  ”;ir.n  2 

99:  wrt  "A.  6",  ',[1,01 

100:  next  T;wth  '1,10,13: next  J 

101:  next  V 

102:  wtb  "4,13 

103:  ate  " 

104;  one 
*31767 
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